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(See glossary 


at the end of 
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ard 21 51 s2 


845.0 


2 trad: 


48.4 48.4 82.0 48.4 48.4 


136.0 
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0.0 
A Weasel 


3.0 
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1 2o 10 
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0.6 
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0.6 
4.7 
156.5 
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4.7 
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1,565.0 


Urea & Carbon 
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145.3 
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14.5 


14.5 


14.5 


14.5 


0.0 
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0.0 
14.5 
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0.0 
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EMISSIONS BY SCENARIO (tonnes/year) 


CONTAMINANT 
(See glossary 


at the end of 
the table.) 


PLANT TYPE 


SIC 


CODE 


A 


Existing 
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9 

6 = 5S 

237.0. 
76 


328.0 


3925 


237..0='2,409.8 
76.0 


237.0 
o.. 


2,412.0 
328.0 


Acid 


1c 


PM 
HC 
wees 
Ni 


Ammonium Nitrate 


(99.5%) 


Anmonium Nitrate 


(63%) 


76.0 


AGT 


109% 


76.0 


712.0 


IM 
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240.6 
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350.0 
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NOx 
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Total 3711 General 
“ Chemicals 
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PM 
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EMISSIONS BY SCENARIO (tonnes/year) 


CONTAMINANT 
(See glossary 


at the end of 
the table.) 


PLANT TYPE 


SIC 
CODE 


A 


Existing 


TOTAL INORGANIC CHEMICAL, 


15.2 


1Se2e5 168-4 
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1552 
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EMISSIONS BY SCENARIO (tonnes/year) 


CONTAMINANT 
(See glossary 


at the end of 
the table.) 


PLANT TYPE 


SIC 


CODE 


A 


Existing 


0.0 
3.6 


24.1 
3369.0 


Ink 


Total 3751 Paint & Varmish 


3791 


a9 1.3. 2.0 


0.1 


1.3 


143 


34.4 


VOC 


Solvent 


0.0 


0.0 0.0 


0.0 


0.2 


Pigment 


0.0 0.0 0.0 0.0 0.4 


0.4 


VOC 


Adhesives 


3792 


Other Chemicals 


Metallic 


SMS 


5.4 
190.0 
498 .3 


5.4 


182.7 
498.3 


35.0 
182./ 
996.7 


5.4 
i2h.3 
498.3 


5.4 


LO 25-7. 
498.3 


52.4 

190.0 
9,700.0 
4,012.0 


Plasticizer & 


19 


Solvent Reclamation 


Pigment Paint 


659.4 


659.4 


659.4 1,310.9 


659.4 


659.4 


659.4, 
182.7 
498.3 


659.4 1,310.9 
12S 
498.3 


659.4 


4,012.0 


190.0 
498.3 


182.7 
996.7 


LOZ 
498.3 


190.0 
9,700.0 


PM 
NOx 


Total 3799 Other Chenicals 


VOC 
Pb 


5.4 5.4 356.0 5.4 5.4 


52.4 


TOTAL ORGANIC CHEMICAL PRODOCTS 


890.9 


890.9 


890.9 3,577.2 


890.9 


6,350.3 


8,639.8 5,102.7 2,846.8 5,462.6 5,102.7 5,114.5 


26 ,;29440-> 2,052.) )2,052.1-6-3 7901.0 2,0 ete oe) 


0.1 
5.4 


0.1 0.1 0.1 
35.0 5.4 


0.1 
5.4 


130755 


5.4 


52.4 


EMISSIONS BY SCENARIO (tonnes/year) 


OONTAMINANT 
(See glossary 


at the end of 
the table.) 


PLANT TYPE 


SIC 
CODE 


A 


POWER GENERATTON 


Electric Power Industries 


Lakeview 


A911 


996.4 5,13024 1,971.6 5,130.4 
841.3 6,247.1 1,664.7. 6,247.1 


7,515.6 3,292.2 SU, 663.8" 07,575.60 G,292. 2" a7, 57566 


5,130.4 1,971.6 
6,247.1 1,664.7 


500.4 
300.0 
443.9 


370.8 
194.0 
443.9 


500.4 
300.0 
443.9 


187.2 


370.8 
194.0 
443.9 


500.4 
300.0 
443.9 


98.0 
443.9 


Raa ee 


Lambton 


Nanticoke 


Thunder Bay 


Atikokan 


Lennox 


= 


TOTAL POWER GENERATTON 


20, 19734. 97, 93722 4, 23036 20,1974 2/93 722 20719754 


20 


WASTE DISPOSAL INDUSTRIES 


0.0 4.5 0.0 0.0 


0.0 


50.0 


Odour Units 


Transfer Station 


4999 


67.3 
0.3 
3756 
0.2 
1.2 
0.3 


67.3 
0.3 
3136 
OZ 
2 
O.3 


101.0 
125 
56.4 
0.8 
1.8 
0.9 


67.3 
0.3 
37.6 
G2 
Loe 
0.3 


67.3 
0.3 
37.6 
OS2 
a De 
0.3 


3,298.0 
34.0 
1,843.0 
19.0 
58.0 
0.9 


Municipal Landfill 


EMISSIONS BY SCENARIO (tonnes/year) 


CONTAMINANT 


PLANT TYPE 


SIC 


CODE 


(See glossary 
at the end of 


the table.) 


A 


Existing 


240.9 240.9 


240.9 
307.0 


S-S Incinerator 


90.4 
180.7 
120,5 


90.4 
105.4 


90.4 
105.4 
120.5 


105.4 
120;5 


105.4 - 
220; 5 


188.0 


120.5 


571.0 


2.4 
0.2 


2.4 
0.2 


264 
Gee 


2.4 2.4 
Go2 


0.2 


2.4 
0.2 


pS DS 2,5 2.5 eS 
0.6 


0.6 


2.5 
0.6 


0.6 


0.6 


0.6 


0.0 
0.0 
0.2 


0.0012 
0.0023 


O22 


0.0012 
0.0023 


0.2 


0.0012 
0.0023 


0.2 


0.0012 


0.0012 


0.00230 0.0023 


0.2 


0,2. 


aul 


Bohs 7 35 Lae 7 35% 7 stag eas 357.7 


SOT T 


MC Incinerator 


25.6 25.6 32.0 2546 25.6 
102.0 


32.0 
102.0 


1647 


1641 


76.7 


16.54 


17.0 17.0 42.6 17.0 17.0 


42.6 


0.01 
0.01 
0.22 
0.01 
0.0 


0.0 


0.03 
0.01 
1.10 
0.06 
0.0 


0.01 
0.01 
0,22 
0.01 
0.0 


0.01 
0.01 
0.22 
0.01 
0.0 


0.03 
0.01 
1.10 
0.06 
0.0 


0.0 
2 


0.0 
0.0 


0.0 0.0 0.0 0.0 
0.8 


0.0 


0.0 


0.06 0.06 0.30 0.06 


0.30 


EMISSIONS BY SCENARIO (tonnes/year) 


CONTAMINANT 
(See glossary 


at the end of 
the table.) 


PLANT TYPE 


SIC 
CODE 


A 


4999 


29.3 29 3 29.3 29 «3 29.3 


15255 


MH Incinerator 


35ee 
176.0 


35 x2 35 a2 3D 28 35 <2 
176.0 B22 oo 176.0 


176.0 


SS ai; 
703.8 


58.25 58.5 58.5 58.5 5855 


53.25 


11.7 Lay. 1M bat Lies 11.7 


17.6 


0.02 0.02 0.02 0.0 


0.02 


0.9 


29353 293 293 2a03 293 


106.7 


FB Incinerator 


28.4 
LOS) 


28.4 28.4 28.4 28.4 
103.5 152.4 105.5 


105.5 


28.4 
152.4 


a2 


58.6 58.6 98.26 58.6 58.6 


58.6 


1). Lia Lal lls Ley 


17.6 


0.02 0.02 0.02 0.02 0.0 


0.10 


Private Waste Disposal 


4999 


4.9 
40.6 


4.9 
40.6 


4.9 
40.6 
1L2heJ 


4.9 
40.6 


4.9 
40.6 


4.9 
40.6 
124.5] 


Industrial Incinerators IM 


Commercial and 


SO2 


Ta 73.0 73.0 
1785 


7320 
118.5 


178.5 


1785 


178.5 


178.5 


24.3 24.3 24.3 24.3 24.3 


24.3 


0.0 
0.2 


0.004 
0.2 


0.040 
0.2 


0.004 
0.42 


0.004 
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0.040 
0.2 


EMISSIONS BY SCENARIO (tonnes/year) 


CONTAMINANT 
(See glossary 


at the end of 
the table.) 


PLANT TYPE 


A 


Institartional Waste 


Di 


4999 


0.5 
0.4 
0.3 
OZ 


0.5 
0.4 
0.3 


1.0 
0.7 
0.6 


0.5 
0.4 
0.3 
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0.5 
0.4 
0.3 
O22 


°108.3 


PM 


Hospitals 


Fics 
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Ac 


rganic 
NOx 


O 


58.0 


O.2 


0.5 


46.6 


Aldehydes 


971.0 657.7 657.7 658.2 657.7 657-7 


466.7 
1,325.8 


2205.1 
536.8 


2200 
536.8 


226.5 
780.2 
186.5 


220.1 
536.8 


220.1 
536.8 
161.0 


Zo 


161.0 


161.0 


161.0 


648.8 


2e5 
0.4 


oan: PS 220 259 
0.4 0.4 0.4 
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6 HS 
0.4 


ed 
0.6 
0.0 
0.0 
0.2 


A AD A es 2.0 dy 
0.6 On 7 
106.1 


2.6 
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0.6 


0.6 


0.0012 

0.002 

0.2 
106.1 


0.0012 
0.002 
0.5 

159.2 


0.0012 

0.002 

0.2 
106.1 


0.0012 

0.002 

O22 
106.1 


0.0012 
0.002 
0.5 


5,199.0 


0.8 
0.4 
0.2 
L118 .5 


0.8 
0.4 
OFZ 


0.8 Sez 
178.5 


0.4 


0.8 
0.4 
0.2 


17855 


53.9 


OF}; 


Ce 
46.6 
18.5 


Organic Aci 


O=5 
178.5 


0.2 
118s 


141.6 141.6 141.6 141.6 141.6 


141.6 
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The emission§ ‘data supplied by Senes*Consultants Ltd: (see 
Appendix D) include two broad categories of contaminants, namely 
particulate matter (PM) and volatile organic compounds (VOC's). 
These broad categories include dozens (PM) or hundreds (VOC) of 
different chemicals. | 


Reduced emissions of particulate matter, regardless of the 
composition, can yield health (respiratory problems, eye 
irritationstetes) and other (selling, «visabllity reduction) 
benefits. But most health and environmental benefits derived 
from reduced emissions of PM's and VOC's result from lower 
concentrations of specific chemicals. Thus, the composition of 
the VOC and PM emissions must be determined so that the changes 
in emissions of specific chemicals can be examined. This was one 


of the major challenges of the project. 
E.2 APPROACH 


Not surprisingly, the composition of the VOC's (PM's) emitted 
varies with the specific industrial process. The emissions from 
a coal-fired boiler, a gas-fired boiler and an oil-fired boiler 
are different. To determine the composition of VOC (PM) 
emissions, it is necessary to begin with the processes used in 


each economic sector. 


E.2.1 “Speciation 


The only source of detailed data on the composition of VOC's and 
PM's is of the United States Environmental Protection Agency's 
(EPA), National Emissions Data System (NEDS). 


As in Canada, U.S. industries are classified by Standard 
Industrial Classification (SIC) codes. Unfortunately, U.S. SIC 


codes differ~from Canadian SIC codes. 


In the NEDS database each U.S. SIC is divided into processes 
which are regarded as sources of emissions and are assigned a 
unique Source. Classification: Cods) (SCe eye To use caer oc 


information on PM and VOC composition we needed to: 


. select the ScC's that correspond to each of the Canadian 
industries or industry segments covered by this study; 


- ensure that the SCC's selected correspond to the processes 
and emissions estimated by Senes; and 


aggregate the SCC profiles of PM's and VOC's into profiles 
corresponding to the Canadian SIC codes covered by this 
study. 


These steps are described below. 
Ess IDENTIFICATION TOR SCCts 


The, first task is to identify the unique, Source Classi ficaticd 
Codes (SCC's) needed for our analysis. The Canadian SIC codes 
for the industries or industry segments covered by the study are 
known. Statistics Canada provides a conversion table that 
relates.each Canadian STC code to the corresponding Ugo 16 
codes. (We used, this conversion table to identify the relevar. 
U.S. SIC codes,’ The number of U.S... SiC codes) identiriacws- 
approximately 50 per cent larger than the corresponding number of 
Canadian SIC codes. 


We then obtained from the National Emissions Data System (NEDS) 
data for the SCC's corresponding to the identified U.S. SIC 
codes. We received. data on approximately 15,000 Scec"s. 


E.4 SCC DATA 


The data obtained for each SCC included: 


the U.S. SIC code that uses the process; 

. the type of process; 

- a description of the process including notes such as the 
fuel used; 
the corresponding speciation profile number which identifies 
the specific chemicals that comprise the VOC (PM) emissions 
for that process: and 
the total VOC (PM) emissions from the given process by all 
plants in the United States in 1985. 


In theory each SCC has a unique chemical breakdown for its voc 
emissions and its PM emissions. In practice, generalized 
profiles are common to many SCC's. For example, the composition 
of the VOC and PM emissions from oil-fired boilers might be 


Similar for Many INCUStrles, (Sls) that, use such soLlléers. 


It is also important to understand that the SCC's, corresponding 
to a given SIC code do not represent a single set of processes 
foundsing all plants:dn that» industry .«o<« Stickingi wath the boiler 
example, a given SIC code may have ScC's that correspond to 
oil-fired, coal-fired, and gas-fired boilers. A given plant 
might. use: only asscingle-?uel for its boiler, so only one’ of the 


sc(Cts”" would-be relevant for its emissions. 


E.5 ELIMINATION OF MINOR SCC'S 


be trie 
As mentioned above, a large number (more than 15,000) processes 
(ScC's) are sources of emissions in the industrial sectors we 
considered. We had to find some way of cutting down this number 
to make the calibration of SCC's with the industry sectors 
covered by the study and the emissions data provided by Senes a 


manageable task. 


We used the following procedure. We eliminated all SCC's whose 
contribution to total PM or total VOC emissions in: the (U7s.) 

SIC in question was less than 5 percent.1t In most cases, we were 
left with a small group of .SCo"s which contributed 75° percenceaen 
more of total PM or VOC emissions of U.S. firms in thateinrdus.. = 


Era a few instances we were left with less than 75 percent of 
total PM or VOC emissions using the "less than 5 percent" rule. 
In those cases we conducted an ordered search for "minor" SCC's 
and. outithem backwinstheuSCG qrouping) for the American S107 
question. We kept adding "minor" SCC's until we had 75 percent 
of total PM and total VOC emissions accounted for. Our ordered 
search involved looking for the most important SCC's of less than 
5 percent, i.e., 4-5 percent, the 3-4 percent, etc. 


The process of eliminating minor SCC's reduced the number from 
ovetm 157000 to ~Lessethan: 1.0002 


- Only if the SCC contributed less than 5 percent of both 
total VOC and PM emissions was it excluded. If the Scc 
contributed 5.1 percent of total VOC emissions=and 
1.0 percent of total PM emissions it was retained. 


sche al aR 


E.6 ELIMINATION OF SUB-SECTORS NOT COVERED 
BY “Leto SruD: 


ar. ea 
Rh) 
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The American SIC codes include sub-sectors that are not found in 
Ontario or were excluded from this study. For example, the 
American SIC code for fish processing includes canneries which 
are not found: in Ontario. Hence, the canneries sub-sector and 
related SCC's were excluded. Similarly, only the coating 
operations of the automotive and automotive parts industries are 
addressed in this study. Hence, SCC's relating to other 
processes in the automotive and automotive parts industries were 
eliminated. 


Few problems were encountered in matching the sector and 
sub-sector definitions to those within the scope of this study. 
Matching the process definitions, with the exception of boilers, 
also did not create any problems. Where questions arose we 
contacted Senes Consultants to ensure that the sub-sector and/or 
process coverage in terms of SCC's matched emission sources they 


Taencitled An-ontar.c. 


E.7 BOILERS 


Boilers are an integral part of many industrial processes. They 
are also used to provide space heating for larger buildings. It 
is not Surprising then to fina chat. a number of the SCC's related 
to boilers. However, Senes advised us that boilers are not 
important sources of emissions in many of the sectors studied and 
hence boiler emissions were not included in their estimated 


emissions for these sectors. 


Bhp ae =) 


Boilers therefore presented three problems for the analysis: 


Whether, for a particular economic sector, the boiler 
aspect, should be included or not. 


Whether boilers used for space heat should be considered. 


If it was decided to include boiler emissions ina 
particular sector, how the speciation of Voc and PM 
emissions obtained from the U.S. should be adjusted to 
account for the greater use of natural gas and oil in 
Ontario relative to the Unitea states. 


Our approach to these problems is described below. The approach 
is apparently simplistic, but we considered it better to be 
consistent in a project with uncertainty, in any case, so 


prevalent. 


E.7.1 Boiler Inclusion in the Analysis 


The EPA NEDS database may indicate that boilers are sources of 
voc's and PM's ina particular U.S. industrial sector. 

Regardless of that information, we followed the advice of the 
consultants who provided the emissions data (Senes) as to whether 
boilers had in fact been included in their calculations. If 
boiler emissions were included in the emissions estimates, EPA 
boiler speciations were used for PM's and VOC's emitted by the 
sector in question. If boiler emissions were not included in the 


emissions estimates, the SCC's relating to boilers were dropped. 


Eliminating the emissions from boilers had the effect of 
reducing the percentage of total VOC or PM emissions covered by 
the selected SCC's below 75 percent of total voc and PM emissions 
by thatoindustrytansmost cases. 


E-7 


Biota Boiler Use for Space Heat 


Emissions from boilers used for space heating were excluded in 
all cases, except to the extent that the boiler would be included 
as part of the general boiler use in the sector. 


E.7.3 Boiler Fuel Conversion, Sectors 


We judged it necessary to make boiler fuel conversions in three 
sectors: electric power generation, pulp and paper, and 
petroleum refineries. 


Ontario Hydro information for the six power generation plants 
examined in this study indicates that only three fuels are used: 
lignite coal, bituminous coal, and heavy fuel oil. The EPA data 
also indicated boilers fueled by wood bark waste, distillate oil 
and natural gas. We treated these fuels as heavy fuel oil and 


adjusted emissions as described below. 


In the pulp and paper sector, we assumed that wood bark waste, 
natural gas, and heavy fuel oil are used in Ontario. Where EPA 
identified the use of coal, we treated these emissions as if 
they came from heavy fuel oil. For petroleum refineries, we 
assumed natural gas and heavy fuel oil are used to fuel boilers. 
Where EPA referred to use of light oil, we converted this to 


heavy fuer ot: 
E.7.4 Boiler Fuel Conversion, Procedure 


Next we had to convert emissions from one fuel type to another. 


We used the following procedure. 


Senes provided us with information.on VOC and PM emissions per 
tonne of steam for different fuels by boiler type (i.e., the many 
different kinds of boiler that are powered by one type of fuel). 


Bree 


voc and PM emissions per tonne of steam do not vary greatly by 
boiler type. The emissions per tonne of steam depend much more 
on the fuel used than on the boiler design. We calculated 
typical voc, and PM emissions per tonne of steam by fuel type, 
regardless of the boiler type. 


To convert the EPA data from one boiler fuel to another, we 
started from the premise that the steam required would not 
change. Then the quantity and composition of the voc and PM 
emissions resulting from the generation of that steam were 
adjusted to reflect the different boiler fuel. 


For example, assume the NEDS data for the pulp and paper industry 
indicated that particulate emissions from bituminous coal-fired 
boilers amounted to 500 tonnes per year. The steam generated was 
estimated at (500 x 1,000/5.35 =) 93,458 tonnes per year.? 
Assuming this same 93,458 tonnes of steam was generated by a 
heavy oil fired boiler, the PM emissions would be (93,458 

Xx 0.075/1,000 =) 7 tonnes per year.? The magnitude of this shift 
in PM emissions -- from 500 tonnes to 7 tonnes in the example -- 
indicates why we felt it important to adjust the mix of boiler 
fuels where possible. 


The PM emissions estimated by converting from one boiler fuel to 
another (bituminous coal to heavy oil) were added to the 
emissions reported by the NEDS database for the substitute fuel 
(heavy oth). s*The «<SCC*profilesfor*boilens using thateinuelwwae 
assumed to apply to the revised total emissions. 


2 500 tonnes of PM x 1,000 kg/tonne divided by PM emissions 
Of 5.35 KO per, tonne Omrevean. 


393,458 tonnes of steam x 0.075 kg of PM emissions per 
tonne of steam/1,000 kg per tonne = 7 tonnes of PM 
emissions. 
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These calculations eliminate some SCC profiles -- those for the 
boiler fuels that have been converted -- from the set of profiles 
for the industry and assign higher total PM and voc emissions to 
other SCC profiles -- those for the substitute fuels ==6. for that 
industry. The total emissions of PM and VOC by SCC obtained from 
the NEDS database, adjusted for the different mix of boiler 
fuels, are used as weights for the SCC profiles in the 
calculation of overall PM and voc emissions profiles for the 
industry. 


This process may appear unsatisfactory. However, the reader 
should remember that the objective of the speciation process is 
to define the characteristic proportions of contaminants in the 
emissions from a particular economic sector. Thus, apparent 
errors resulting from our adjustments may alter the composition 
of the VOC or PM emissions speciations somewhat, but will not 
affect the total amount of such emissions. The amounts are those 


estimated by Senes. 
E.8 CALCULATION OF SPECIATIONS 


Having» ddentitiedithessScetsithat aecountitoretne bulkeoé The ®Vvoc 
and PM emissions in the sectors and sub-sectors of interest for 
this study and having adjusted the process fuels where 
appropriate, it remained to calculate a VOC and PM speciation 
for each industry. That was done as follows. 


All of the selected SCC's corresponding to a given Canadian SIC 
were identified. The composition of each of these SCC's was 
determined from the NEDS database. This lists every chemical 
that is found in the VOC (PM) emissions by that Scc together with 
the percent weight of that chemical of total VOC (PM) emitted by 
that processwSCC)S by. UsSclplantsoin 19959 aNexcy weighted 


EO 


averages: of the same contaminant from different SCC's were 
combined. This total was used to determine each contaminant's 
percentage share of the total voc (PM) emissions for a specific 


Canadian Si CIR « 


In a few instances, the emissions data provided by Senes 
Consultants (Appendix D) identifies separately chemicals that 
are constituents of the PM or VOC emissions. We checked with 
them to determine whether those chemicals were emitted 
separately, in addition to the emissions that of that chemical as 
part of the PM or voc emissions. Where the emissions were said 
to be part of the PM or VOC emissions, the quantities of that 
chemical estimated by the speciation of the PM or VOC emissions 
were subsequently replaced by Senes' estimate of emissions for 
the contaminant. This uses the more relevant information 
generated by Senes, and eliminates double counting of emissions 


Of those’ contaminants. 
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The PM and VOC speciations for the Canadian SIC codes used in 
this analysis are shown below. 


The profiles break down the PM and voc emissions into their 
constituent chemicals. Some of those constituent chemicals are 
Regulation 308 Contaminants. Others’ are not on the listvcr 

96 chemicals we were asked to consider. Each profile shows the 
percentage contribution of each Regulation 308 Contaminant and a 


balance percentage of "other" chemicals making’ up: the» profive: 


This "other" category typically accounts for most -- usually over 
60 percent -- ofthe PMser VOCy emissions’) inital sectormmanmre 


relatively large size of this "other" category raised the concern 
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that potentially harmful contaminants were being overlooked. To 
allay this concern we provided a list of all. voc and PM chemicals 
on the NEDS database to the Ministry of the Environment for 
review. The Ministry was satisfied that the list ofpeontaminancs 
originally specified identified the chemicals of concern. This 
review led to the addition of three chemicals to the list of 


contaminants. 


Thus, the "other" category can be considered to consist of 
relatively harmless chemicals. The steps taken to reduce the 
more hazardous chemicals in the VOC or PM emission streams will 
have the effect of lowering emissions of these more benign 
constituents as well at no additional cost. Because of the 
relatively benign nature of these "other" chemicals no benefits 
are deemed to be derived from their curtailment. 


pM BOP 
OPEN PIT MINING 
SICOGGold Mines 
PM Contaminant “8 Der Cent 


SILICA 31.42 


OTHER 68.58 
Subtotal 100.00 


SIC 0617 Iron Ore Mines 


IM Contaminant Per Cent 
ARSENIC 2.64 
CADMIUM 0.75 
CHROMIUM Oe2s 
FLUORIDE G52 
MERCURY 0502 
MANGANESE 0.68 
NICKEL 0.14 
LEAD 4.93 
SELENIUM 0.05 
SILICA 4,54 
AGIOS | Oe20 
OTHER So.oU 
Subtotal 99.99 


FOOD INDUSTRIES 


SIC 1011 Meat ard Meat Proecucts 


EM Contaminant Per Cent VOC Contaminant Per Cent 
SILICA 15.00 BENZENE Osti 
OTHER 85.00 ae: ETHYLBENZENE 0.01 
Subtotal 100.00 ETHYLENE . 0.63 
. FORMALDEHYDE 0.81 
TOLUENE 0.04 
OTHER 98.40 
Subtotal 100.00 


SIC 31012 ) Pouliry Products 


EM Contaminant Per Cent . voc Contaminant Per Cent 
SILICA 15.00 OTHER 100.00 
OTHER . 85.00 Subtotal 100.00 


Subtotal 100.00 


Eo-i13 
SIC 1021 Fish Preducts 


FM Contaminant Per Cent 


SILICA Sa ae 

OTHER scooeacia8 57-00 

Subtotal 100.00 

SIC 1049 Other Dairy Products (Whey 

Drying) 

FM Contaminant Per Cent voc Contaminant 
ARSENIC 03°25 OTHER 
CADMIUM O07 Subtotal 
CHROMIUM 0.03 

FIUORIDE 0.05 

MERCURY 0.01 

MANGANESE 0.07 

NICKEL 0.03 

LEAD On6Sm 

SELENIUM 0.01 

SILICA lsiw6 

TIN ereQ3.03 

OTHER 85.04 

Subtotal LOG nO 

PULP AND PAPER AND ALLIED PRODUCTS 

SIC 7271... Pulp: Indus=ary 

PM Contaminant Per Cent Voc Contaminant 
ARSENIC O02: - BENZENE 
CADMIUM 03.0% ETHYLBENZENE 
CHROMIUM OF 0S : ETHYLENE 
FLIDORIDE 1.19 FORMALDEHYDE 
MERCURY “ (OF.0S TOLUENE 
MANGANESE On.15 OTHER 
NICKEL 07,03 Subtotal 
LEAD 01.02 

SELENIUM 05,02 

SELICA OF 55 

TIN 0.02 

OTHER 97.90 


Subtotal 100.00 


pol ag 


100. 
100. 


é 
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Cent 


00 
00 
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SIC 2712 Newsprint Industry 


PM Contaminant Per Cent Voc Contaminant Per Cent 
ARSENIC Pes cca OPA 1,1,1-TRICHLOROETHANE 0.02 
CADMIUM Sine ROO 1, 3-BUTADIENE 0.06 
CHROMIUM 0.02 ACRYLONITRILE 0.05 
FLUORIDE 0.85 BENZENE 9.46 
MERCURY 0.05 BUTYL ACRYLATE 0.03 
MANGANESE O12 BUTYL CELLOSOLVE C502 
NICKEL 0.43 CARBON TETRACHLORIDE 0.04 
LEAD 0.02 CELLOSOLVE 0.02 
SELENIUM 0.01 CELIOSOLVE ACETATE 0.02 
SILICA 203 CHLOROFORM 0.03 
TIN 0.02 EPICHLOROHYDRIN 0.02 
OTHER 96.43 ETHYLBENZENE O73 
Subtotal 100.00 ETHYLENE 7.48 
ETHYLENE DIBROMIDE 0.02 
ETHYLENE DICHLORIDE "0.09 
ETHYLENE OXIDE 0.02 
FORMALDEHYDE 4.55 
METHYL ETHYL KETONE 0.11 
METHYL ISOBUTYL KETONE 0.04 
NAPTHALENE 0.03 
PERCHLOROETHYLENE 0.04 
PHENOL 0.55 
PROPIONALDEHYDE 0402 
PROPIONIC ACID 0.02 
PROPYLENE OXIDE 0.03 
STYRENE Onl? 
TOLUENE 2 ee 
TOLUENE DIISOCYANATE 0.04 
TRICHLOROETHYIENE 0.02 
VINYL CHLORIDE 0.06 
OTHER 74.06 


Subtotal 100.01 


1S ogee =) 


SIC 2713 Paperboard Industry 
FM Contaminant Per Cent 


ont 


ARSENIC ' eieyor 
CADMIUM pa oe eet COROT 
CHROMIUM 0.01 
FIUORIDE 1.18 
MERCURY 0.08 
MANGANESE 6.17 
NICKEL 0.50 
LEAD | 0.02 
SILICA 2.06 
TIN 0.02 
OTHER 95.92 
Subtotal 99.98 


SIC 2811 Business Fomns Printing 


SIC 2819 Other Comercial Printing 


EM Contaminant Per Cent 


ARSENIC 1.04 
CADMIUM 0.34 
CHROMIUM 0.12 
FIUORIDE 0.22 
MERCURY 0.05 
MANGANESE 0.32 
NICKEL 0.12 
LEAD 3.10 
SELENIUM 0.03 
SILICA 8.36 
TIN O15 
OTHER 86.17 


Subtotal 100.02 
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SIC 2821 Plate Making, Typesetting ard 


Bindery 

FM Contaminant . | Per Cent 
CADMIUM , GC. O7 
FIUORIDE 0.36 
MANGANESE 0.04 
NICKEL GO. Ou 
LEAD O30L 
SILICA Oe2n 
TIN Oe 
OTHER 99.335 
Subtotal 100.00 


SIC 2841 Newspaper, Magazine and 


Pericdical Publishing 


ETHYLENE OXIDE 
FORMALDEHYDE 

METHYL ETHYL KETONE 
METHYL ISOBUTYL KETONE 
NAPTHALENE 
PERCHLOROETHYLENE 
PHENOL 
PROPIONALDEHYDE 
PROPIONIC ACID 
PROPYLENE OXIDE 
STYRENE 

TOLUENE 

TOLUENE DIISOCYANATE 
TRICHLOROETHYLENE 
VINYL CHLORIDE 


. OTHER 


Subtotal 


i 
e 
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OOO0 00 00 00 020 00 00 0 0 RnA0N0000 


Ear Lis. 


SIC 2849 Other Combined Publishing and 
Printing Industries 
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IRON AND STEEL INDUSTRIES 
SIC 2911 Ferro-Alloys 


FM Contaminant Per Cent 


ARSENIC 0.21 
CADMIUM 0.20 
CHROMIUM Lace 
MANGANESE 0.29 
NICKEL 0.10 
LEAD 3.56 
SELENIUM 0.02 
SILICA 6.61 
TIN 0.01 
OTHER 87.61 
Subtotal hese be 


Subtotal 


ETHYLENE DICHLORIDE 
ETHYLENE OXIDE 
FORMALDEHYDE 
METHYL ETHYL KETONE 
METHYL ISOBUTYL KETONE 
NAPTHALENE 
PERCHLOROETHYLENE 
PHENOL 
PROPIONALDEHYDE 
PROPIONIC ACID 
PROPYLENE OXIDE 


é 


}? 
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100. 


Cent 


- 00 
-00 
-00 
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SIC 2912 Steel Foundries 


PM Contaminant § Per Cent 
ARSENIC a LO OZ 
CADMIUM i at BO aCe 
CHROMIUM 236 
MANGANESE 3555 
NICKEL 1.04 
LEAD 0.50 
SILICA 2.92 
TIN 0.01 
OTHER 90.75 
Subtotal 99.98 


Voc Contaminant 


1,1) i-fRIGALORORTHANE 
1,3-BUTADIENE 
ACRYLONTTRILE 
BENZENE 

BUTYL ACRYLATE 

BUTYL CELLOSOLVE 
CARBON TETRACHLORIDE 
CELLOSOLVE 
CELLOSOLVE ACETATE 
CHLOROFORM 
EPICHLOROHYDRIN 
ETHYLBENZENE 
ETHYLENE 

ETHYLENE DIBROMIDE 
ETHYLENE DICHLORIDE 
ETHYLENE OXIDE 
FORMALDEHYDE 

METHYL ETHYL KETONE 
METHYL ISOBUTYL KETONE 
NAPTHALENE 
PERCHLORO: ; 
PHENOL 
PROPICONALDEHYDE 
PROPIONIC ACID 


Per Cent 


0.43 
0.67 
0.58 
2.82 
0.43 
0.44 
0.65 
0.36 
0.36 
0.51 
0.42 
0.65 
2.86 
0.42 
0.58 
0.38 
1.55 
hele 29 
0.56 
0.18 
0.70 
0.44 
0.42 
0.40 
0.50 
»i20 
2.04 
0.66 
0.01 
0.43 
0.42 
76.67 
99.96 


SIC 2919 Other Primary Steel Industries 


FM Contaminant 


al a 


CADMIUM 
CHROMIUM 
FIUVORIDE 
MANGANESE 
NICKEL 
LEAD 
SELENIUM 


- Per Cent 


\ 
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100.00 


FOUNDRIES 


SiG) 2940. Iron Foundries . 


MM Contaminant 


ARSENIC 
CADMIUM 
CHROMIUM 
MANGANESE 
NICKEL 
LEAD 
SILICA 
OTHER 
Subtotal 


“i\ Per Cent 


Com 
O07 
0.08 
3.04 
0.10 
0.38 
13.58 
82.59 
99.99 


voc Contaminant 


1, pt TRICE ORORTHANE 
1,3-BUTADIENE 
ACRYLONITRILE 
BENZENE 

BUTYL ACRYLATE 
BUTYL CELLOSOLVE 
CARPON TETRACHLORIDE 
CELLOSOLVE 
CELLOSOLVE ACETATE 
CHLOROFORM 
EPICHLOROHYDRIN 
ETHYLBENZENE 
ETHYLENE 

ETHYLENE DIBROMIDE 
ETHYLENE DICHLORIDE 
ETHYLENE OXIDE 
FORMALDEHYDE 

METHYL EPHYL). KETONE 
METHYL ISOBUTYL KETONE 
NAPTHALENE 
PERCHLOROETHYLENE 
PHENOL 
PROPIONALDEHYDE 
PROPICNIC, ACID 


TOLUENE. 

TOLUENE DIISOCYANATE 
TRICHLOROETHYLENE 
VINYL CHLORIDE 
OTHER 

Subtotal 
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SIC 2971 Copper amd Copper Alloy 
Rolling, Casting and Extruding 


FM Contaminant = Per Cent 
ARSENIC  OSO3 
CADMIUM 0.01 
CHROMIUM 0.02 
MANGANESE alvAl 
NICKEL 0.02 
LEAD 0.37 
SILICA 3.60 
TIN 0.01 
OTHER 95.25 
Subtotal 100.02 


SIC 2999 Other Rolled, Cast amd Extruded 
Non—Ferrous Metal Products 


EM Contaminant Per Cent 
ARSENIC 0.06 
CADMIUM 0205 
CHROMIUM 1.35 

_ MANGANESE OroL 
NICKEL O55 
LEAD 0.95 
SILICA 4.80 
TIN 0.06 
OTHER 91.08 


Subtotal eho ep, 


PERCHLOROETHYLENE 
Subtotal 


00 
00 


eae 
NON-FERROUS METAL SMELTING AND REFINING 


SIC 2959 Other Primary Smelting ard 
Refining of Non-Ferrous Metal 


PM Contaminant - ‘Per cént voc Contaminant Per Cent 
ARSENIC 25.48 BENZENE 4.01 
CADMIUM 1.46 _ ETHYLENE O2OL 
CHROMIUM OvZ2 FORMALDEHYDE 7.98 
FLUORIDE Ls4y7 TOLUENE 2.3.00 
MERCURY 0.04 OTHER 85.99 
MANGANESE Oia Subtotal 99.99 
NICKEL 0.09 

LEAD 1225 

SELENIUM 0.23 

SILICA . 2225 

TIN Oss2 

OTHER 55.66 

Subtotal 99.99 


MOTOR VEHICLE INDUSTRIES 


SIC 3231 Motor Vehicles amd Motor 
Vehicle Parts amd Accessories Surface 


Coating 

FM Contaminant Per Cent VOC Contaminant Per Cent 
CADMIUM 0.2.05 BENZENE 1.04 
CHROMIUM 0.05 ETHYLBENZENE 8.68 
NICKEL 0265 TOLUENE 20 ji2 
OTHER 99.85 OTHER 70.2% 


Subtotal 100.00 Subtotal 100.00 
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SIC 3254 Motor Vehicle Steering ard 
Suspension Parts 


FM Contaminant °. Per Cent voc Contaminant Per Cent 

OTHER Sy geass 1,1, 1-TRICHLOROETHANE 1.89 

Subtotal 99.99 1,3-BUTADIENE 0.96 
ACRYLONITRILE 0.01 
BENZENE 0.39 
BUTYL CELLOSOLVE 0.76 
ETHYLBENZENE 23 
ETHYIENE © 2229 
ETHYIENE DICHLORIDE 0.02 
FORMALDEHYDE 0.44 
METHYL ETHYL KETONE O217 
METHYL ISOBUTYL KETONE 0.04 
PERCHLOROETHYILENE 0.63 
PHENOL 0.01 
PROPIONALDEHYDE 0.01 
STYRENE 0.04 
TOLUENE 22.76 
TRICHLOROETHYTENE 1.79 
VINYL CHLORIDE 0.01 
OTHER 75.45 
Subtotal 100.00 


INORGANIC CHEMICAL MANUFACTURING 


SIC 3511 Domestic Clay Preducts 


EM Contaminant Per Cent Voc Contaminant Per Cent 
ARSENIC 0701 BENZENE ST lo) 
CADMIUM 0.01 ETHYLBENZENE 3229 
CHROMIUM eel ETHYLENE 0.08 
FLIUORIDE O03 FORMALDEHYDE 2.96 
MANGANESE On.07 NAPTHALENE O2xEL0 
NICKEL 0.06 TOLUENE aie, 
LEAD Ges OTHER 84.21 
SELENIUM 0.06 Subtotal 99.99 
SILICA dito l ; 

OTHER 91.76 


Subtotal 100.01 


SIC 3512 Imported Clay Products 


IM Contaminant 


- Per Cent 
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0.04 
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SIC 3521 Hydraulic Cement 


FM Contaminant ‘Per Cent yor Contaminant 

CADMIUM ees OL 1,1, 1-TRICHLOROETHANE 

CHROMIUM LePrecatos 1; 3-BUTADIENE 

FIUORIDE 0.12 ACRYLONITRILE 

MANGANESE 0.02 BENZENE 

NICKEL 0.03 BUTYL ACRYLATE 

LEAD - 0.83 BUTYL CELLOSOLVE. 

SELENIUM 0.01 CARBON TETRACHLORIDE 

SILICA 5.48 CELLOSOLVE 

OTHER 93.29 CELLOSOLVE ACETATE 

Subtotal 100.01 CHLOROFORM 
EPICHILOROHYDRIN 
ETHYLBENZENE 
ETHYLENE 


ETHYLENE DIBROMIDE 

. ETHYLENE DICHLORIDE 
ETHYLENE OXIDE 
FORMALDEHYDE 
METHYL ETHYL KETONE 
METHYL ISOBUTYL KETONE 
NAPTHALENE 
PERCHLORO 
PHENOL 
PROPIONALDEHYDE 
PROPIONIC ACID 
PROPYLENE OXIDE 
STYRENE 
TOLUENE 
TOLUENE DIISOCYANATE 
TRICHLOROBENZENES 
TRICHLOROETHYLENE 
VINYL CHLORIDE 
OTHER 
Subtotal 


f 
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3571 Abrasives 
EM Contaminant 


ARSENIC 


e 


yee So 
FPOUkNOUUN 
; ct 


OFOOCO0000 


Ey="26 
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Sic 3581"" Lame 
FM Contaminant 


ARSENIC 
CADMIUM 
CHROMIUM 
FIUORIDE 
MANGANESE 
NICKEL 
LEAD 
SELENIUM 
STLICA 
OTHER 
Subtotal 


SIC 3711 Industrial Inorganic chemicals 


(NEC) 
EM Contaminant 


‘Per Cent 


eos OF 


Onl. 
0.17 
0.03 
0.07 
0.06 
0.17 
0.06 
LeOL 

Po ae fo) 

100.01 


Per Cent 


100.01 


oe? 


PROPIONALDEHYDE 
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SIC 3721 Chemical Fertilizer ard 
Fertilizer Materials 


FM Contaminant | ‘Per Cent 


« 
Eves 


a 


ARSENIC aie at ORC oe 
CADMIUM 0.05 
CHROMIUM 0.03 
FIUORIDE 0.01 
MANGANESE 0.09 
NICKEL 0.04 
LEAD 0.16 
SELENIUM 0.03 
SILICA 4.04 
TIN 0.02 
OTHER 95.46 
Subtotal 997.95 


PROPIONALDEHYDE 


Per Cent 


0.01 
0.02 
36.80 
2.73 
0.. 07 
0.01 
0.10 
0.01 
0.01 
0.03 
0.01 
0.19 
Sanu 
0.01 
0.69 
0.03 
0.78 
0.55 
0.08 
0.01 
0.02 
0.45 
0.01 
OF02 
0.02 
25.44 
0.73 
0.02 
0.01 
0.44 
0.01 
28.35 
997795 


E - 29 
PETROLEUM REFINERTES 
SIC 3611 Petroleum Refineries 


FM Contaminant “Per Cent 


ARSENIC 0.01 
CHROMIUM 0.14 
FIDORIDE 0.02 
MANGANESE 0.03 
NICKEL Lod 
LEAD 0.07 
SELENIUM 0.13 
SILICA 8.57 
OTHER 89.5L 
Subtotal EE), 


ASPHALT PAVING INDUSTRY 


SIC 3699 Asphalt Paving Industry 
FM Contaminant Per Cent 


ARSENIC 0.01 
CADMIUM 0.01 
CHROMIUM 0.15 
FLUORIDE 0.02 
MANGANESE 0.07 
NICKEL 0.05 
LEAD 0.15 
SELENIUM 0.06 
STLICA 6.73 
OTHER 92975 
Subtotal 100.00 
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ORGANIC CHEMICAL PRODUCTS 


SIC 3712 Industrial Organic Chemicals 


FM Contaminant: - Per Cent 
ARSENIC OX 
CADMIUM 0. 
CHROMIUM 0. 
MERCURY 0% 
MANGANESE 0% 
NICKEL ‘om 
LEAD 0. 
SELENIUM ok 
SILICA 4. 
TIN a 
OTHER 94, 
Subtotal 100. 


SIC 3731 Plastics and Synthetics Resins 


EM Contaminant Per Cent 


OTHER 93. 
Subtotal 100. 
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02 
45 
05 
84 
Ol 


TOLUENE DIISOCYANATE 
TRICHLOROETHYLENE 
VINYL CHLORIDE 
XYLENE BASE ACIDS 
OTHER 

Subtotal 
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SIC 3751 Paint and Varnish 


FM Contaminant | ‘Per Cent 


ek ee 


CACMTIUM 07-03% 
CHROMIUM 0.02 
MANGANESE On 
NICKEL O02 
LEAD Os.O2 
SELENIUM 0702 
SILICA S693 
TIN 0.03 
OTHER 93.81 
Subtotal 99.98 
SIC 3791 Printing Ink 
FM Contaminant Per Cent | 
CALDMIUM 0.04 
CHROMIUM 0.02 
MANGANESE 0.12 
NICKEL 0.03 
LEAD 0.02 
SELENIUM 0.02 
SILICA 5.95 
TIN 0.06 
OTHER 93.80 


Per Cent 


10.45 
0.01 
0.01 
0.01 
hee Od. 
52.56 
6.97 

31.70 

37.84 

99596 


Per Cent 


3.00 
5.00 
92.00 
100.00 
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SIC 3792 Adhesives 


FM Contaminant ‘Per Cent 
CADMIUM | Olas 
CHROMIUM ne eee Oa 
MANGANESE 0.04 
NICKEL 0.04 
SELENIUM Or On) 
SILICA 1.98 
TIN O.02 
OTHER 97.84 
Subtotal 100.01 


SIC 3799 Other Chemical Products (NEC) 


PM Contaminant Per Cent 
CADMIUM 0503 
CHROMIUM OOS 
MANGANESE 0. 0S 
NICKEL Oy 
LEAD OO 
SELENIUM He) 
SLEICA O05 
TIN OF 02 
OTHER 96.70 
Subtotal 100.02 


Voc Contaminant 


1,3-BUTADIENE 
ACRYLONITRILE 
BENZENE 

BUTYL ACRYLATE 
BUTYL CELLOSOLVE 
CARBON TETRACHLORIDE 
CELLOSOLVE 
CELLOSOLVE ACETATE 
ETHYLBENZENE 
ETHYLENE 

ETHYLENE DIBROMIDE 
ETHYLENE DICHLORIDE 
ETHYLENE OXIDE 
FORMALDEHYDE 
METHYL ETHYL KETONE 
METHYL ISOBUTYL KETONE 
NAPTHALENE 

PHENOL 

STYRENE 

TOLUENE 

VINYL CHLORIDE 
OTHER 

Subtotal 


Voc Contaminant 


ACRYLONITRILE 
BENZENE 

BUTYL ACRYLATE 
CARBON TETRACHLORIDE 
ETHYLBENZENE 


.ETHYLENE 


ETHYLENE DICHLORIDE 
ETHYLENE OXIDE 
FORMALDEHYDE 
METHYL ETHYL KETONE 
METHYL ISOBUTYL KETONE 
PHENOL 

PROPYLENE OXIDE 
STYRENE 

TOLUENE 

VINYL CHLORIDE 
XYLENE BASE ACIDS 
OTHER 

Subtotal 
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ELECTRIC POWER SYSTEMS 


SIC 4911 Electric ‘Power Systems - Coal 
and Heavy Oil Fired Generating Stations 


FM Contaminant Per Cent Voc Contaminant 


ARSENIC 0.05 1,3-BUTADIENE 
CHROMIUM 0.06 BENZENE 
MERCURY 0.01 ETHYLBENZENE 
MANGANESE 0.04 ETHYLENE 
NICKEL O25 FORMALDEHYDE 
LEAD 0.17 PROPIONALDEHYDE 
SELENIUM 0.02 TOLUENE 
SILICA 20.06 OTHER 

TIN 0.01 Subtotal 
OTHER ~ Toes a | 

Subtotal 99.29 


WASTE DISPOSAL INDUSTRIES 
SIC 4999 Hospital Incinerators 


FM Contaminant Per Cent 


ARSENIC 0.01 
CADMIUM 0.06 
CHROMIUM 0.02 
MERCURY 2.67 
MANGANESE 0.04 
NICKEL 0.01 
LEAD ; A Slee 
TIN 0.44 
OTHER 92.64 


Subtotal 100.01 


Per Cent 


0.28 
0.37 
5.88 
1.87 
14.19 
oPueRs 
Pdi hs) 
74.39 
Sp hee hs) 
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SIC 4999 Mmicipal and Commerical 


Incirerators 

FM Contaminant . Per Cent Voc Contaminant Per Cent 
ARSENIC Hi iene aes BENZENE 7.28 
CADMIUM 0.09 ETHYLENE Seay 
CHROMIUM 0.03. OTHER 83.43 
FIUORIDE 0.02 Subtotal 99.98 
MERCURY Sie 

MANGANESE 0.07 

NICKEL 0.01 

LEAD 6.15 

SILICA 0.05 

TIN 0.66 

OTHER 88.93 

Subtotal 100.02 


SIC 4999 Sewage Sludge Incinerators 


EM Contaminant Per Cent voc Contaminant Per Cent 
CATMIUM O5S3 BENZENE 74g0 
CHROMIUM 0.58 ETHYLENE 83.70 
MANGANESE Gad OTHER | 83:60 
LEAD 2.94 Subtotal 100.00 
SELENIUM OLS 

SILICA pte AP 

TIN roo 

OTHER 89.13 

Subtotal 100.00 


DRY CLEANING ESTABLISHMENTS 


SIC 9721 Dry Cleaners 


EM Contaminant Per Cent voc Contaminant Per Cent 
ARSENIC be ws) BENZENE O03 
CADMIUM 0.38 ETHYLBENZENE Ope 
CHROMIUM eas FORMALDEHYDE 0307 
FLUORIDE 0.24 NAPTHALENE Oc ie 
MERCURY O05 PERCHLOROETHYLENE BBs ee 
MANGANESE 0.34 PHENOL 4.93 
NICKEL Onn Ss TOLUENE Or Ze 
LEAD 3000 OTHER 81.20 
SELENIUM O03 Subtotal 100.03 
SILICA 8.67 

TIN Oe 
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he 
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Subtotal ° 100.00 


Seas 
vie 
et ‘ 


APPENDIX F 


Dispersion Modelling 
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Fel, LNORODUCTION 


This appendix describes the dispersion modelling work performed 
to estimate grids for annual average concentration of air 


emissions in‘various parts of Ontario. 


FIs lpmoOpganizdtionyof the work 


The principal steps in the work are: 
. selection of a model; 
selection of locations for which the model is run; 
- compilation of required meteorological data for each 
ROCaLI On, 
. specification of source characteristics; and 
- quality assurance checks of the model results. 


The approach taken to each of these steps is outlined in the 
following sections. 
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Lapfbradizza- 


We also appreciate the assistance of staff of the Canadian 
Climate Centre, Atmospheric Environment Service, Environment 


Canada-.andeUsS fhNationab’Chimatic Data ‘centre: 
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F.2 SELECTION OF A DISPERSION MODEL 


rey a Proposéd Dispersion Model 


Our original proposal described a regulatory dispersion model 
presently in use in the Province of Alberta. This model was a 
variant of the existing Ontario regulatory model, (a Pasquill and 
Gifford type). The multi-source, receptor version of the Alberta 
model was considered superior in many respects to the existing 
Ontario regulatory model and thus to be preferred. 


Due to demand generated by widespread use of Pasquill-Gifford 
models, Environment Canada has developed the STAR (STability 
ARray) program to generate the simplified climatological input 
data used by these models. Typical input consists of frequency 
distributions of wind direction, wind speeds, and the Pasquill 
stabilityilc lasses (se Vs). Gy iete 


F.2.2.° Ministry of the Environment's Models 


In discussions with the Ontario Ministry of the Environment's 
Atmospheric Model Development Unit, during the early stages of 
this project, it was agreed that use of models currently under 
development in Ontario would be preferred. The models under 
development were far superior to the Pasquill-type models. The 
only relative shortcomings of the newer models were in the ease 
of use. For example, meteorological input: data are in’ the form 
of detailed hourly meteorological observations. This represents 
a quantum leap in the quantity of data to be handled and gave 
rise to considerable complications in data transfer and data 
management. 


A development version of the Ministry's multi-source, (GAS) model 
was released to the study team in early June. This helped 


ae 3 


facilitate data conversion program development. A final 
"pre-public-release" version was released to the study team on 
June 21, 1988.'. That version was used to execute the modelling 
runs. eres 


Some minor modification to the meteorological parameter 
calculation program (MET) provided by the Ministry was undertaken 
to flag additional meteorological cases as invalid (so these 
cases would not be calculated). Program and data conversion 


development are described in more detail in subsequent sections. 


A full discussion of the new models is presented in Appendix H of 
the Clean Air Program Discussion Paper issued by the Ontario 
Ministry of the Environment in late 1987 (Environment Ontario, 
1987). 


F.3 SELECTION OF LOCATIONS 


Fe3cees Tntrodtiiction 


To determine the public benefits of the proposed revisions to 
Regulation 308,.it is necessary to determine the dispersion of 
emissions of each of the 96 specified contaminants emitted by 
each of the over 3,500 sources identified. In principle, each 
source must be modelled separately because the source 
characteristics and meteorological conditions are unique. [In 
practice that is impossible due to lack of data and the magnitude 
of the effort it would require. 


The availability of climatological information allows only a 
limited number of sites to be modelled; those for which the 
required input data can be compiled. The dispersion patterns 


modelled for those locations are then applied to all sources 


Peeig 


within the area to which those climatological characteristics 
apply. This, in turn, requires that standard characteristics he 
assumed for the emission sources. The use of standard source 
characteristics is not a major concern since the required 
information is, with very rare exceptions, not available in any 
case. The source characteristic assumptions are discussed in 
Section F.5 below. Here we review the procedures used to 
identify the locations with different climatological 


characteristics. 


F.3.2 Climatological Data Requirements 


The Ministry's GAS model requires input climatological data from 
both hourly observing stations and upper air stations. There are 
no locations within the study area for which both hourly 
observations and upper air data are available. To get the model 
input required upper data from U.S. and Quebec upper air stations 
must be matched with data from Ontario hourly observing stations. 


The process of selecting locations to be modelled starts with the 
list of hourly observing stations in Ontario and eliminates those 
that have similar climatological conditions. Appropriate upper 


air data are then matched to these remaining locations. 


The Ministry's MET and GAS models require hourly air temperature, 
cloud opacity, wind speed and wind direction observations in 
addition to a daily snow on the ground indicator and mixing 
height, which is derived from upper air temperature profiles. 
Thus ‘the starting point for thevselection, of locations sis then2 
hourly observing stations in Ontario in the climatological 
archive of Environment Canada, Atmospheric Environment Service. 
Most of these stations have a data record extending for 
approximately 30 years. This is the typical period of record 
considered necessary for a useful climatological data set. 


baie 12) 


One year is composed of 8760 hours, so the analysis of four 
hourly weather elements for a 30 Year period wfor all “Ontario 
stations would require analysis of over 30 million observations. 
This was well “beyond the scope of the Dro Veci. 


P.3e3 Approach to the: SéLection sotsLocations 


To reduce the number of locations, and hence the computational 
work load, in a scientifically credible manner it is necessary 
to develop criteria to determine when stations have similar 
Climatological characteristics from the perspective of the 
dispersion model. 


An examination of the Ontario Ministry of the Environment's model 
indicated that the development of dispersion contours around a 
source depends heavily on the wind speed and direction of the 
source area. Accordingly, we began by comparing mean annual 
statistics of wind speed and direction for the 32 available 
stations. 


Several stations are located very close to one another, and so 
allow us to define criteria for similar climatological 
characteristics. For example, Toronto Downsview Airport 
(latitude 43° 45' N, longitude 79° 29' W) or Toronto 
International Airport (latitude 43° 40' N, longitude 79° 38' W) 
could be used to represent conditions in the Toronto area. 
Similarly, either Ottawa Rockcliffe Airport (latitude 45° 27', 
longitude 75° 38' W) or Ottawa International Airport (latitude 
45° 19' N, longitude 75° 40' W) could be used to represent 
conditionsrin the Ottawa arean™ Ic follows that’ if other 
geographically proximate stations had wind statistics as similar 
to one another as these station pairs, then they too could be 


considered interchangeable for the. purposes of this study. 


Er eS 


F.3.4 Wind Speed and Direction Comparisons 


Stations were first compared by examining the average annual 
frequency of “otcurrence of each wind direction (8 point compass) 
and calm conditions. Each of these nine values for one station 
in a station pair was subtracted from the corresponding value for 
the other station to obtain their difference. The largest 
difference across these nine values was then selected to 
represent the maximum departure of each station from the other in 
terms of wind direction frequency. For the two Toronto stations, 
the largest difference occurred for northerly winds which are 
observed 17 per cent of the time at the International Airport as 
opposed to 11.7 per cent of the time at the Downsview Airport 

for a difference of 5.3 per cent (Environment Canada, 1982). The 
corresponding value for the two Ottawa stations was 

3.4 per, cent. These values are given tin eexhibit oF akvalong mien 


Similar values for all possible station pairs in Ontario. 


Wind speed characteristics were compared by examining average 
annual wind speed values by direction (8 point compass). Again 
the largest difference among the eight wind speeds was selected 
to represent the maximum departure of each station from the other 
for each station pair. These values are given in Exhibit F.2 for 
all Ontario station pairs. . \Roriche two Toronto, stats cons ibe 
maximum wind speed difference was 2.7 km/h, while the 
corresponding value for the Ottawa stations was 3.4 km/h. 


FPvgw5y Sebvecttonvcrivteria 


From these examples we concluded that stations could be 
considered to be essentially similar for dispersion modelling 
purposes if their maximum wind direction difference was less than 
5 per cent and if their maximum wind speed difference was less 
than 3° km/heceThisyrule is, ifianvthing: ss tooestringenmes nce 
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Eee US 
neither the Toronto stations nor the Ottawa stations would 
qualify as essentially similar according to our criteria. 
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Foo «6 Select'ién efuLocations 


Four geographical areas were found to contain "essentially 
Similar" stations that met or almost met the similarity criteria 
defined above. These areas included northwestern Ontario, 
northeastern Ontario, Georgian Bay and southwestern Ontario. 
Plots of average annual wind direction occurrence and average 
annual wind speed by direction for the candidate stations in 
these regions are given these regions in Exhibit F.3 through F.6, 
respectively. The two stations in Georgian Bay and Southwestern 
Ontario obeyed the similarity rule completely. 


In northeastern Ontario, Timmins Airport and Kapuskasing Airport 
were essentially similar according to the rule while the Timmins 
Ailrport/Earlton ALrport pair (Saceper centrteant4pe-—imy ny and 
the Earlton Airport/Kapuskasing Airport pair (3.8 per cent and 
3.21 km/h) did not. In northwestern Ontario, Armstrong Airport 
was the station most closely approximating the four others in 
this group,..but,.onLy.the- Armstrong Airport/7Ssioux Lookout pair 


were essential similar according to our selection rule. 


Despite the fact that some northeastern stations and most 
northwestern stations had differences larger than was specified 
by the criteria of essential similarity, these stations are 
located in regions with very few establishments in the economic 
sectors being studied. Hence the use of additional dispersion 
patterns in those regions was not justified, despite the slightly 
larger dissimilarities. 


The amalgamations adopted based on the procedures described above 
are as follows: (the representative station is underlined) ; 


Toronto (Toronto International Airport, Toronto Downsview 
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Wind Correlations-Northwestern Ontario 
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Wind Correlations-Northeastern Ontario 
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Wind Correlations-Georgian Bay 
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Wind Correlations-Southwestern Ontario 
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Airport), Ottawa (Ottawa International Airport ;-0ttawa-Rotkelitie 
Airport), Northwestern Ontario (Armstrong Airport, Atikokan, 
Geraldton, Pickle Lake, Sioux Lookout), Northeastern Ontario 
(Timmins Airport, Earlton Airport, Kapuskasing Airport), 

Georgian Bay (Wiarton Airport, Gore Bay Airport), and 
Southwestern Ontario (London Airport, Mount Forest). These 
amalgamations resulted in the removal of 10 stations from the 


data base. 


Eight additional stations were also removed since they are not 
in areas with significant sources of the contaminants studied in 
this project. These stations are Armstrong Airport, Big Trout 
Lake, Lansdowne House, Moosonee, Muskoka Airport, Petawawa 
Airport, Red Lake Airport, and Wawa Airport. Any sources 
examined in these regions use more distant stations to define 


average meteorological conditions. 


The stations retained, which together are taken to define the 
meteorological conditions for the province, are Kenora Airport, 
Kingston Airport, London Airport, North Bay Airport, Ottawa 
International Airport, Sault Ste. Marie Airport, Simcoe, Sudbury 
Airport, Thunder Bay Airport, Timmins Airport, Toronto 
International Airport, Trenton Airport, Wiarton Airport and 
Windsor Airport. 


F.4 COMPILATION OF METEOROLOGICAL DATA 


F.4.1 Need to Select a Representative Year 


While the rationalization of observing stations to be included in 
the analysis reduced the data to be analyzed by more than 50 per 
cent, the data burden was still too large for the project. To 
accomplish a further reduction it was decided to replace the 
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30 year period of record used to characterize meteorological 
conditions with a representative year that closely approximates 
the 30 year period. This representative year yields an annual 
dispersion fiéld very close to that obtained from the 30 year 
history for each of the 14 stations. 


For the purposes of the project only annual statistics 
summarizing the dispersion pattern are required, the large number 
of hourly observations in one year ensure that statistics derived 
from these hourly observations do not vary greatly from year to 
year. Consequently, any year in the data base could be used to 
represent typical conditions if all input data were hourly 
observations. Such, however, is not the case. 


F.4.2 Snow#cover 


Specifically the presence or absence of snow on the ground is 
not subject to the law of large numbers in the way that hourly 
observations are. The period of time that snow is present for 
any station varies a great deal from year to year. In the 
Ministry's model the presence of snow on the ground is assumed 
Cor retrect ait~sticort=wave radiation lpack™= into space; = wich, che 
result that convective conditions cannot occur (Ontario Ministry 
of the Environment, 1987). The disallowance of convective 
conditions for an unusually long or unusually short period 
relative to the 30 year average may well result ina 
Significantly different dispersion field. “Accordingly, it was 
decided to seek a single year which, for all 14 stations, the 
winter snow cover period was as close as possible to the average 


obtained using the entire period of record for the station. 


The time period for which snow was on the ground for at least 
seven consecutive days was determined for each of the 14 
stations. for all years of each stations's climatological record. 


Significant year to year variability was found. For example, in 
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Kenora a 30 year data base indicated that on average seven 
consecutive days of snow cover began on November 11 and ended on 
Aprerl 11 for . a ‘total of 152 days. However, in 1968 snow was 
present for. only 126 days while in 1965 snow was present for 179 
days. Other stations showed similar variability in terms of the 


length of the snow cover season. 


F.4.3 Selection of the Representative Year 


Only one year, 1986, had the time with snow on the ground within 
one standard deviation of the long term average for all 

14 stations. Thus 1986 was selected as the representative year 
for the calculation of average annual dispersion statistics. 


F.4.4 Compilation of Climatological Input Data: Hourly Station 
Data 


Hourly meteorological observations for the representative year 
for each of the 14 centres (totalling 500,000 observations) were 
purchased from the Canadian Climate Centre, Atmospheric 
Environment Service and formatted for input to the Ministry's 
METCO program. 


The METCO program, provided by the Ministry, allows the use of 
data formatted for use in the single source GAS model in the 
newer multi-source, multi-receptor version of the model. The 
program also allows disabling of upper air wind speeds where 
desirable. 


This program was used as the basis for a general conversion 
program (METCOR) written by the study group to automate the 
interface of Environment Canada Surface Data, U.S. and Canadian 
Upper Air Data, the snow record, data decision-making, flagging 
and reformatting of input data. The METCOR program produced 
14 complete records for the representative year, ready for input 
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to the Ministry's MET program (meteorological parameter 
Calculation program). These 14 records represent average 
meteorological’ conditions fOr, tne Province: or Ontariovror the 
period 1957, to'1986 inclusive. The MET program generates the 
meteorological parameters required as input by the dispersion 
(GAS) model. 


F.4.5 GMT Conversion 


Canadian surface weather data are recorded in local standard time 
for each station. A program segment was written to convert local 
standard time (date) to Greenwich Mean Time (GMT) for input to 
the GAS model. 


Fae ormcloudecover 


Canadian surface weather data expresses cloud cover and cloud 
opacity in 1l0ths. Cloud cover and cloud opacity affect the 
occurrence of convective conditions. Cloud opacity, where 
available, is the better determinator of convective conditions 
and was used in this study. A program segment was used to 
convert from 10ths to 8ths for input to the GAS model. 
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The most significant determinators of convective conditions are 
the presence or absence of snow cover and the presence or 
absence of daylight. Most winters are characterized by 
alternating periods of "snow-on" and "snow-off". Wide swings 
occur, throughout the 30-year record, of starting and ending 
dates - dates of snow cover, and dates of bare ground. The snow 


cover period may also be interrupted by periods of bare ground. 


The 30-year record for each centre was reduced to a mean start 


and vend*“date- of*snow cover.  “The’def—inition of snow cover used is 


YR Rasa =: 


the period from which persistent (7-day duration), 2 cm snow 
cover was first recorded to its last recorded date. This 
definition had the effect of reducing the length of the snow 
season. The’ “Feduction in the length of the season offsets, to 
some degree, the increase in the number of "snow-on" days which 


result from calling the snow cover period "continuous". 


The 30-year mean period was used as the "snow-on" period for the 
representative year. A program segment was written to automate 
snow cover distribution in the model. 


Canadian surface wind speed is in units of km/hr. A program 
segment was written to convert this to m/sec for the GAS model. 
Calm conditions were assigned a wind speed value of 0 m/sec where 
no value was provided. The MET and GAS modules increase this to 
a 1.5 m/sec minimum in appropriate cases. Upper air wind speeds 
were disabled in the model. 


FP.4.9 Wind Direction 


The representative year chosen, being recent, was available in 

10 degree compass resolution and was thus compatible with the 

GAS model interpolation routine (wind direction +5 degree 
randomization). A program segment was written to multiply the 
AES compass, (0 to.36) by 10.. In,Environment Canada format...3621s 
North and 0 is calm. A program segment was written to reassign 
36 to 0, for._the GAS model. Initial, Oo (calm) values. were 

assigned the first previous non-calm wind direction. 


F.4.10 Upper Air Data 


No Canadian upper air stations exist in proximity to most 
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Ontario emission centres. Thus, it was necessary to obtain U.S. 


near-border data to use with Ontario surface data. 


Pred ols) Mixing‘Heights: UsiowrLocatirons 


Mixing heights were purchased for the representative year from 
the’ U.S. National Climatic Data Centre’ in North’ Carolina. © The 
U.S. computer program used, generates a.m. and p.m. mixing 
heights using upper air soundings combined with surface data for 
the same station (Doty, 1983). The normal U.S. practice is to 
combine surface data from the surface location of interest. 
However, since the surface locations of interest are all 
Canadian centres, this was not possible. The use of matching 
upper air/surface data for the same station was considered a 


reasonable compromise. 


Twice daily, mixing heights for the entire year were obtained 
for, Flint, Michigan; «Sault, Ste. MariefeMtchacgans~ international 
Falls, Minnesota; and Buffalo, New York. 


Fi4 l2i Mixing Heights canadtan tocations 


Since mixing heights cannot be purchased from Environment Canada 
for Canadian Upper Air Stations, raw soundings data were 
purchased. The only Canadian centre in proximity to the 
locations to be modelled is Maniwaki, Quebec. 


Since Maniwaki, Quebec surface data were not purchased, Ottawa 
surface data were used with the Maniwaki upper air data to 
generate mixing heights. This is in accordance with U.S. 
practice. A program was written to adjust base temperatures for 
base elevations, determine maximum surface temperatures for each 
period, and calculate the elevation of intersection between the 
sounding temperature profile and the dry adiabatic lapse curve 


for each surface temperature maximum. 


Fe ae 


A linear interpolation routine was written to refine 
intersections which occurred between sounding elevations. 
Negative elevations. were assigned a value of zero, which was 
flagged by the revised MET program for cases which would 
ordinarily be convective. In such cases, a zero mixing height 
resulted in an invalid MET case in the GAS model execution 


(ves; nercalecularionin, 


Ottawa was chosen as the surface centre because it is 
geographically central to the area of Ontario for which 
Maniwaki, Quebec upper air data would be deemed representative 
(i.e., Timmins, North Bay, Ottawa, Kingston). Time and budget 
did not permit a centre-by-centre match of local surface and 


Maniwaki upper air data. 


F.4.13 Missing Data Flags 


In the case of missing meteorological data a flag value of IDATA 
= 999 was assigned. This generates an invalid MET case for the 


GAS run“ (isespnoecalculation for that case). 


F.4.14 Negative or Missing Mixing Height Flags 


In cases of negative or missing mixing heights, IDATA was 
assigned a value of 987. A short addition to the MET program was 
written such that if IDATA = 987 and meteorological conditions 
were convective, IFLAG would be assigned a value of 999 (an | 
invalid case). 


F.4.15 Limitations of the Meteorological Database 


Some centres, despite catalogue indications, were found to have 
incomplete meteorological records. For example, in Simcoe, 


measurements were only available for weekday periods. Results of 
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modelling with such data are less reliable than results for more 
detailed input. 


The use of Upper Air Data from geographically removed centres is 
a somewhat dubious procedure. However it is one that is without 
alternative since upper air stations are so rare. For example, 
using Buffalo upper air data for Toronto leads to questions of 
how representative upper air data for a centre on the east end of 
a lake is of another centre on the north shore of a different 
lake. Errors would occur only in the case of convective 
conditions and are estimated to yield an output variance, on an 
annual average, of no more than 5 to 10 per cent. The best 
centre, by this criterion, is Sault Ste. Marie, Ontario, since 
the Sault Ste. Marie, Michigan upper air station is so close to 
the Canadian surface weather station. | 


There are many surface stations in northwestern Ontario. Kenora 
was chosen to represent this area based on its population 
density. It is estimated that this choice would affect the 
dispersion of pollutants in other centres (e.g., Armstrong) to a 
smal ler degree (Srtosie” per cent).-"it°- was not expected” that 
large concentration excursions would occur and none were 
observed. 


The dispersion model (GAS) does not deal with calm conditions "as 
are" but rather increases wind speed to 1.5 m/sec. for applicable 
cases. In addition, Environment Canada merely flags wind 
direction in calm cases (since by definition of calm, there is no 
wind to have a direction). To mesh what the model does, with 
what the source data doesn't provide, a wind direction was 
assigned in each case of calm. The choice of direction was that 
of the first previous "non-calm" wind direction. This perturbs 


the short term wind direction Statiscrc Yess “than” i1ntrodtetion of 
an.artificial randomizgationefunction &,oine,gréequency er 


occurrence of calm conditions is estimated at less than 
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2 per cent, so the impact of this assumption on results 


considered to be small. 
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F.5 EMISSION SOURCE CHARACTERISTICS 


The GAS model requires, in addition to climatological data, input 
with respect to emission source characteristics. The 
characteristics of the several thousand emissions sources are, 
with very rare exceptions, not known. Hence these 
characteristics must be estimated. Since each of the 

14 dispersion patterns could be applied to several hundred 
sources, the estimated characteristics should be typical of the 


variety of industrial establishments emitting the contaminants. 


F°5.1. Building Sources 


Of the several thousand commercial emissions sources in Ontario 
subject to the proposed revisions to Regulation 308, by far the 
most numerous are building sources. Such sources are generally 
defined as sources from which fugitive emissions emanate or stack 
emissions emerge at elevations less than twice the height of the 
bud lean. 


To precisely model emissions from such sources the exact 
dimensions of the building and the stack must be known. In cases 
of complex building shape (i.e., non-rectangular) or complex 
vicinity (Lie.,..largeribuildings..or variablesterrain) 
mathematical modelling is inappropriate and physical modelling is 
preferred. Physical modelling is not possible for this analysis. 
Hence it is necessary to assume that all buildings are isolated, 
rectangular buildings with a stack less than twice the height of 
the building. 


1D eae 5) 


For this study, no information was available on building or 
building-stack dimensions for any source. To estimate ambient 
concentration,’ each source was assumed to be a building of 7 m 
height and 100'm width with a stack of less than 14 m elevation, 
emitting pollutants at 0 velocity and at ambient exit 
Penpermacure. such a HUuLlUIng 1s cd, typical, “sndustridl pan” 
source. This assumption leads to very high ambient 
concentration uncertainties (i.e., order of magnitude) within 
500 m of the source for source buildings of other configuration 
but is a far better source assumption than a zero elevation 
emitter. Reasonable concentration estimates beyond 500 m are 
achieved. 


Since no data were available on individual elevated and buoyant 
(exit temperature above the ambient temperature) elevated source 
characteristics, such sources have been represented as building 
sources. The limitations of this approach are discussed in the 


next section: 


F.5.2 Implications of the Assumed Building Source 


Characteristics 


If an emission is ground-based, the building source assumption 
wilicivd eldvatrtut rcPapliyenicneanitra, drinktomwothin oOCrm cfiiithe 
origin. If an emission is elevated, and possibly buoyant, the 
building source assumption will yield artificially high 
concentrations within 500 m of the source. If the source is less 
than about 30 m in elevation and of moderate buoyancy, beyond 

500 m concentrations will generally fall within the same order of 
magnitude with the highest discrepancies found at the perimeter 
(near the zero concentration area). In general, the higher and 


hotter the source, the larger the discrepancies. 


Distribution inaccuracies resulting from the assumptions can be 


illustrated by the fact that the worst case meteorology for an 
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elevated source occurs under convective atmospheric conditions, 
whereas the worst case meteorology for a building source occurs 
under stable. atmospheric conditions. This type of possible 
inaccuracy: céutributed to the decision to use concentration 
averages of 8,760 hourly meteorological cases rather than 
averages over shorter time frames. Where sources approached the 
characteristics: of the "typical. building, source" the: accuracyvior 
the annual average would be unassailable. In other cases, the 
annual average would have a "smoothing" effect where periods of 
larger differences would be offset to some degree by periods of 
smaller differences yielding a better estimate of average 
concentration. This minimizes the limitations of these source 


assumptions. 
F.5.3 Elevated Sources 


Elevated sources, by definition, emit beyond twice the elevation 
of the highest part of an attached or adjacent building. 
Depending on the stack exit temperature the emission plume may be 
more or less buoyant than the surrounding air. Accurate 
modelling requires that all of the stack parameters be known. 
These parameters are: stack height, exit diameter, exit 


velocity, exit temperature (and emission rate). 


Such information was not collected for any source in Ontario as 
pact iL whhiis study. In addition, the use the proposed new 
Ontario models entirely precluded use of individual source 
information of any kind. If one took 2,000 sources for example, 
and ran one years! meteorology, using an IBM-AT with 
math-coprocessor, 36,000 hours of computer time would have been 
required (i.e., 4.1 years). 
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F.5.4 Near-Lake Sources 


Discussions with Ministry isteff led to consideration that 
lakeshore effects should be dealt with in modelling certain 
sources. 


Predominant sources to which lakeshore effects would most often 
be encountered were identified as Lakeview Generating Station in 
Mississauga, and Nanticoke Generating Station. The present 
development version of the Ministry's proposed Lakeshore Model 
has not been integrated with the general GAS model. Such 
integration is necessary because a lakeshore effect can only 
occur under specific conditions which will only arise at certain 
times during the year. At other times dispersion would occur as 
per the GAS model. : 


We endeavoured to perform the model integration; however, the 
project budget was exceeded as a result and this work remains 
incomplete. Hence, lakeshore effects are not included in the 


oubputwort thisstudys 


The error involved in ignoring lakeshore effects for these two 
sources would be an under-estimation in the range of 1 to 3 times 
of annual average concentrations within 3 km onshore from these 


two sources. 
F.6 THE RECEPTOR GRID 


The Ministry's grid-generator program RECGEN was used to generate 
a receptor grid network from 0 to 24 km of 2 km resolution. A 


finer grid of 500 m and 1,000 m grid points was superimposed on 
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the coarse grid to permit better resolution in the near-source . 


area where most of the impact of building sources is felt. 


pers 
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The GAS model ‘calculated running average (annual average in this 
study) concentrations for each receptor point (grid coordinate). 
A surface roughness of 1.0 m and a receptor elevation of 0.0 m 


was used in all cases. 


Due to the assumptions relating to source characteristics the 
resulting ambient concentrations within 500 m are most 

uncertain. As an approximation, where receptor grids 
intersect/overlap (for different sources),if the source is known 
to be elevated the within-500 m area could be taken to be of zero 
concentration. Where the source was not known to be elevated, 
the within-500 m area could be taken to be of concentration equal 
to that at the 500 m boundary. The alternative would be to 
ignore all areas of overlap within 500 m of the source (which in 
the case of Toronto might mean ignoring most of the city for 

some pollutants). 


Fourteen dispersion patterns were generated; one for each 
location selected. The source characteristics are assumed to be 
the same for all 14 patterns. An emission rate of 1 gram/sec is 
used in each case. Since the model is linear with respect to 
dispersion, the calculated concentration at each grid coordinate 
is scaled by the emission rate for the specific source. For 
example, if the concentration at twelve kilometers due east of 
the source is calculated by the model to be 0.0158 micrograms/ 
cu.m. and a given source emits a contaminant at a rate of 

5 grams/sec then the ambient concentration twelve km due east of 


that usource .will 4bé@ 0:07.00 microcrams Accum. 
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F.7 QUALITY ASSURANCE 


F.7.1 InputQuality Control 


Considerable processing, conversion, and flagging of data were 
required prior to running the dispersion model. Conversion, 
processing, and flagging programs were all validated manually 
using Toronto and Simcoe meteorological data. MET program 
output (GAS program input) was reviewed in detail (i.e., for the 
entire one year hourly record) prior to GAS program execution. 
Byte volumes of records for other centres were compared for 
consistency and random hourly records were inspected for other 
centres. | 


All centres used the same input files, except for meteorological 
data where individual centre files were used. To maintain file 
control, alk non=.MET files? (1Re.FeeRUNG .FUN, BREC, SRC and 
blank .AVG files) were located in a write-protected directory. A 
DOS batch-program was written which transferred these files 
(renamed) to the execution directory, prompted for the correct 
~MEDerile Pwerifived. the) presence: off fcheseist illest mand) executed the 
GAS model program. i 


Fi. eet Out par Quality Contro} 


The number of valid and invalid (uncalculated) MET cases was 
inspected on completion of each run and compared with the prior 
knowledge of the relative occurrence of missing and flagging 
Gatar.o Exiriletey 7, the Runimecornd” ” iigstrateslithis:; 


All output (.AVG). files were imported to a LOTUS Symphony file, 
compressed, and compared. The coordinates of the maximum 


concentration receptor were inspected in each case. 
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The relative concentrations of the maxima were compared. With 
each centre using the same building source and unit emission 
rate, the maxima should be within the same order of magnitude. 

cn addition the concentrations at every grid point were summed 
for each centre and the sums compared between centres. This test 
on distribution should also yield results in the same order of 
magnitude. 


All of the above tests were passed. Exhibit F.8 presents the 
first page of this summary file as an example. 


A DOS batch program was written to copy the appropriate files to 
diskette for back-up and to save output files to the appropriate 
directory. The DOS programs ensured that file integrity was 
maintained during movement and handling. 
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G.1. INTRODUCTION 


A vital Gio] gin the analysis of the public benefits of proposed 
revisions to“ Regulation 308 used was DPA's spatial analysis 
system (SPANS) . SPANS is a geographic information system. 
Different geographic information systems specialise in performing 
different tasks; some work best at creating cartographic maps. 
SPANS' specialty is analysis and manipulation of spatial data; 
data which relate to points or areas on maps. It was ideally 
suited to this assignment which involved locating sources of 
emissions, the dispersion of those emissions, and identification 


of the population and environmental resources affected. 


This appendix describes some of the features of SPANS and the way 


it was used in the project. 
G.2 WHAT IS SPATIAL ANALYSIS? 


Spatial analysis is the ability to investigate the spatial 
dimension of a prablem. Many problems have a spatial dimension. 
For instance, the analysis of the public benefits of reduced air 
emissions depend upon where contaminants are deposited and the 
concentration at each point because that determines the resources 
affected. Traditionally, the spatial dimensions of problems are 
addressed by making assumptions (e.g., trade area boundaries) or 
using manual procedures (map overlays) that limit the number of 


options considered. 


DPA's computerized spatial analysis system allows the spatial 
components of a problem to be addressed in much more 
sophisticated and’ convenient ways. To understand how the 
spatial analysis system works Sone background -informatroin-on 


geographic information systems is useful. 


tert GEOGRAPHIC INFORMATION SYSTEMS 
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Computerized geographic information systems (GIS) are designed to 
input, store, retrieve, manipulate and display spatial 
information. They differ primarily in the manner in which they 
store data and the kind of manipulations they perform best. [In 


turn these capabilities determine the size of computer needed. 


G.3.1l: Data Storace 


A spatial feature has a location, size and shape. The size and 
shape can range from a point (e.g., a meterological station) to 
an irregular area (e.g.., a watershed). in addition to these 
spatial characteristics, the GIS must store a range of attributes 
associated with the site: average annual temperature, soil type, 
slope, population, suitability for swimming, etc. Clearly, the 
data storage demands can quickly become very large, creating a 


need for larger computers. 


The most common approach to accurately specify spatial features 
is to..define their. boundaries... A,series of line segments are 
linked to form the boundary of an area. Obviously, the more 
irregular the area the larger the storage requirements. Some 
analytical processes, such as map overlays, are very complex 
with these systems. Line segments must be used to generate the 
areas, the areas are overlaid, the portions that overlap are 
determined and a new set of boundary line segments is generated 


for storage. These systems typically require mini-computers. 


The usual approach to minimising storage needs is to pre-define a 
fixed set of areas, Typically :thissiswa.grid (raster es] ome 


although other systems, such as census subdivisions, can be used. 


GS 


The fixed areas and locations defined by such systems yield major 
benefits in terms of data storage. The cost is loss of accuracy. 
Watersheds, slopes, soils, rainfall and many other attributes do 


not vary according to a fixed grid pattern. 


SPANS uses a quadtree structure to store data. This offers a 
combination of the best features of a boundary-based system and 
grid system. The quadtree structure divides the total area into 
progressively smaller quadrants, until a quadrant has uniform 
charactcristics...Tnis results dn arearunits ot varvabiec Sizes 'so 
a large area with uniform attributes is stored as such, while 
much smaller area units are used to accurately define the 
boundaries of spatial features, such as forest cover, lakes, and 
Dusit. ip area. ghe overall data storage requirements are 
Similar to those for a fixed area system, but the level of 
accuracy is much greater. The system operates on a micro 


computer. 


G.3.2e,ManipulationsGapabilities 


The manipulation capabilities of a GIS system reflect the needs 
Of users and »jtossome sextent jpotherndatatstorage (techniques 
employed. The available systems can be broadly defined as 
cartographic or thematic in orientation. The principal uses of 
cartographic systems are to record, process revisions to, and 
prepare updated maps. Although these systems have some 
analytical capabilities, such as the ability to perform overlays, 
the emphasis is on cartographic properties such as accuracy, 


labelling, changes to scale and projection, etc. 


Conversely thematic systems emphasize analytical capabilities but 
at the cost of some cartographic ability such as labelling. 
SPANS was specifically designed as an analytical system. Its 


@apabriities include: 


Giana 


Easy access to a variety of input data; 
Extensive data manipulation capabilities; 
Complex overlays of multiple maps; 
Proximity, mapping; 

Area analysis reports; 

Use of satellite imagery; and 


Diverse data display capabilities. 


SPANS grew out of needs encountered in the course of DPA's 
consulting work over the past decade. The system was developed, 
refined and enhanced over many years and a wide variety of 


applications. 
G.4. USE, OF .SRANS) IN «THIS sPROJECT 


A critical requirement in the calculation of the public benefits 
of reduced air emissions is the ability to estimate the change 

in concentration of a given contaminant at any point and to 
determine the environmental resources that will benefit from this 


reduced exposure. SPANS provides this | Capability. 
G.4.1 Location of Emission Sources 


SPANS is,.used to locate each emission,source., The list. of 
establishments was compiled from the Made in Ontario data base 
maintained by the Ministry of Industry, Trade and Technology and 
Statistics Canada's list of establishments. The Made in Ontario 
data base includes the postal code and these establishments were 
located, at..the (six digit) postal code centroid. The centroid 
locations for postal codes were purchased from Statistics 


Canada. The Statistics Canada list only gives the municipality. 
These establishments were located at the centroid for the 


ir oro 


municipal unit. The municipality centroid was calculated using 
SPANS. 

plete) facilitate the calculation of depositions from different 
sources, a standard 4 km grid was created. Each establishment 
was moved to the nearest grid point, a distance of less than 

2 km. Given that the initial location was the postal code or 
MUNLCIDaLLCVOCeruL Oi, “Clie moro RiOLe li  LOGuCe way a lLonie Leaner as 
Or error into the analysis. (Grid points with numerous ‘(over 10) 
establishments were examined. The larger firms assigned to that 
grid point were contacted and, if appropriate, shifted to other 
grid points. However, it remains true that establishments are 
artificially concentrated on some grid points due to the lack of 
accuracy of the address information. The establishments from the 
Statistics Canada list are concentrated at fewer locations 


because there are fewer municipalities than postal codes. 


The attributes stored for each establishment are its: 
Canadian SIC code; 
estimated employment (or other size characteristic) ; 
estimated current emissions by contaminant; 
estimated emissions by contaminant for each implementation 
scenario; and 


an indicator to identify the appropriate dispersion pattern. 
G.4.2 Determination of Concentrations 


The dispersion modelling (see Appendix F) calculated 14 
dispersion patterns for different parts of the province. The 
dispersion model assumed an emission rate of 1 gm/sec and 
calculated the resulting concentration (micrograms/cu.m) ona 
2 km grid covering an area 24 km in all directions: from the 


source. 
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To keep the calculations manageable -- 65 contaminants, 

over 3,500 establishments, 5 scenarios -- it was necessary to 
adopt a coarser grid. Thus, ai4 km grid rather than a’2 km grid 
was used. -Siice the sources had been moved to points on a 4 km 
grid, the concentration estimates also corresponded to points on 


this grid. 


The analysis proceeds on a contaminant by contaminant basis. If 
an establishment is a source of the contaminant, the emissions 
for the specified scenario are converted to an emission rate. 

The appropriate dispersion pattern is selected. The origin of 
the dispersion pattern is adjusted to the source location. Every 
point in the dispersion grid is similarly adjusted to points on 
the province-wide grid. The concentration at each point on the 
dispersion grid is scaled by the emission rate and added to the 


existing concentration at that point on the province-wide grid. 


When this process has been completed the grid points reflect the 
estimated concentration of the contaminant at that point due to 
emissions from all establishments under the given scenario. The 
concentration at each grid point is assumed to apply to a 4 km 


square area with the grid point at its centre. 
Gs4.3  Galculation of Ambient weoncentration 


The dispersion model estimates the concentration of a 
contaminant due to a specific source. Many of the contaminants 
analysed are also emitted by sources other than those considered 
in this study. Thus, the ambient concentration due to all 
sources is likely to be higher than the estimated emissions from 


the Regulation 308 sources covered by this study. 


Data on total ambient concentrations are collected by monitoring 
stations. Two networks of monitoring Stations exist in Oncarr. 


One network has stations located in areas where air quality 


Goa! 


problems are believed to exist. The other network has stations 
located more uniformly across the province to collect data to 
monitor acidic deposition. Each monitoring station collects 
information ‘on, ambient concentrations of a number of 
contaminants. Of the 96 contaminants, monitoring data are 
available for cadmium, chromium, lead, manganese, nickel, 
nitrogen oxides, nitrous oxide, sulphur dioxide, total 
particulate matter (PM) andwtotalevocts. 


For each of these contaminants except total VOC's, the ambient 
concentration as recorded by the monitoring station was compared 
with the concentration at the monitoring station location 
estimated using the dispersion model. The result is a series of 
"residuals" (one for each monitoring station) that represent the 
ambient concentration due to non-Regulation 308 sources. These 
residuals were examined individually for each contaminant. Based 
on this examination, one (or more) values of the ambient 
concentration due to non-Regulation 308 sources were established. 
The value(s) and the rationale for the selection is presented in 
Chapter 4 in the section relating to the specific monitored 


contaminant. 


The value(s) of the concentration due to non-Regulation 308 
sources is added to the estimated concentration due to 
Regulation 308 sources at each grid point for all scenarios. 
This converts the value at each grid point from an estimate of 
the concentration due to the specific Regulation 308 sources 
analysed to an estimate of total ambient concentration. And it 
permits an analysis of the impact of the proposed changes to 
Regulation 308 on total ambient concentration of the given 
contaminant. The ambient concentration at each grid point is 
assumed to apply to a 4 km square area with the grid point at 


rts boentre: 


Gre 


This procedure is only possible for the nine contaminants ae 
cadmium, chromium, lead, manganese, nickel, nitrogen oxides, 
nitrous oxide, sulphur dioxide and PM -- for which monitoring 
data are ated tele The procedure was not used WLERSVOC Ls 
because the monitoring stations report total voc's and we did not 
have a reasonable basis for speciating this total at each 


monitoring station. 
G.4.4 Maps of Environmental Resources at Risk 


Having determined the geographic dispersion of a contaminant and 
its concentration at each point, it is necessary to relate that 
information to the resources at risk. These include humans, 
crops, livestock, forests and surface water. Maps .of.these 
resources were purchased from various sources and converted to 
SPANS format. 


Population by enumeration area for the 1986 census was purchased 
from Statistics Canada. The enumeration area population totals 
were aggregated to the census subdivision level. The area of 
each census subdivision was calculated by SPANS. This allowed 
the population figure to be converted toa population density 
(people per sq. km). The population density is assumed to be 


uniform across: the census subdivision. 


Agricultural data from the 1986 Census of Agriculture were 
purchased from Statistics Canada. The information is available 
at the consolidated census sub-division level. Only the total 


crop acreage by consolidated census sub-division was mapped. 


Forestry data were obtained from the National Forestry Service in 
Petawawa. A map of surface water for Ontario was purchased from 


the Canada Lands Directorate Service. 
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G.4.5 Analysis of Environmental Resources at Risk 


BOr most contaminants, only the estimated emissions due to 
Regulation’ 308..sources are available. .In these cases where the 
analysis can only identify the environmental resources at risk. 
The dispersion pattern for the contaminant is superimposed on the 
appropriate environmental resource(s), usually population. This 
is done for each scenario. The difference between the current 
resource exposure and the exposure under each scenario is 
calculated. This indicates the maximum change in environmental 


resources at risk by scenario. 


In each of these cases, the estimated emissions due to 
Regulation 308 sources are also compared to the threshold values 
above which observable effects occur. In virtually every 
instance the current emissions are well below the threshold 


level. 


G.4.6 Analysis of Contaminants for which Ambient Concentrations 
have been Estimated 


The nine contaminants -- nickel, chromium, lead, manganese, 
nickel, nitrogen oxides, nitrous oxide, sulphur dioxide and PM 
-- for which total ambient concentrations are estimated are 


analysed in a manner similar to that outlined above. 


The areas where the total ambient concentration exceeds the 
threshold value for observable effects is plotted for the 
existing situation and for each scenario for implementation of 
proposed changes to Regulation 308. The difference in the 
resources (population, crops, forests) exposed to above-threshold 
levels of the contaminant is a measure of the public benefits of 


the proposed revisions to Regulation 308. 
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G.4.7,. Analysis ‘of ‘Cases with Exposuréeskesponse and Damage 


Functions 


In the case of cancer, and systemic health, visibility and 
soiling effects in the case of a few important contaminants, 
exposure-response and damage functions are available. These 
functions calculate the benefits of reduced exposure from the 


change in exposure and the population exposed. 


For each scenario, the reduction in exposure from the existing 
level is calculated at each grid point. This change is 
multiplied by the population density in the 4 km square of which 
the grid point is the centre and by the area (16 km“) to get the 
population and change in exposure. These values are summed over 
all relevant grid points and the aggregate values are used in the 
exposure-response and damage function to calculate the estimated 


benefits of reduced emissions. 


Maps showing the dispersion patterns by contaminant are provided 
-in Appendix I. 
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H.1. «INTRODUCTION 


An exposure-response function predicts how health and other 
BaACtors, suchas visibility and materials damage, change as the 
concentration of a contaminant changes.  Exposure-response 
functions are the most commonly available for deaths by cancer. 
Functions are also available for other health effects from 
nitrogen dioxide, sulphur dioxide, and suspended particular 
matter. 


The exposure-response functions are discussed below for cancer 
and various effects of nitrogen dioxide, sulphur dioxide, and 
suspended particulate matter. The functions used are listed 
under the contaminant in question in Appendix Cc. They are 
elaborated here for greater clarity. For many non-cancer 
effects, upper, lower, and central values for the function are 
listed in Appendix C. We have calculated the lower and central 


Value for each function. 
jhe CANCER FUNCTIONS 


che change in Canter dearns cue to a change sf the conesntration 
SL a Contaminan® to preatcted, USing’ a uniterisk, factor... The 
unit risk factors we used are usually IRIS (Integrated Risk 
Information System) values. Where IRIS values are not 


available, we used Massachusetts unit risk factors. 


The change in annual cancer mortality from a change in the 
Concentration @f aA chemical rs Calcwiratea using ene following 


formula: 


Hi => 


ACM) = URE 3x6 (CONC 2 POP 740 


Where: 

ACM = Change in annual cancer mortality resulting from the 
eae in concentration of the contaminant. 

URF = The IRIS (or Massachusetts) Unit Risk Factor. 

CONC = Change in concentration of the contaminant measured on 
an annual basis expressed in micrograms per cubic 
metre. 

POP = The population subject to the change in contaminant 
concentration. 

70 = A factor converting lifetime (estimated to be 70 


years) deaths to annual deaths. 


The change in annual cancer mortality (ACM) is calculated for 
each’ contaminant) for: which: ar cancerfunit. risk) factor, is 
available. The calculation is performed for each of the 


implementation scenarios. 


The calculation is performed as follows for a given contaminant 
and scenario. Each of the over 15,000 grid cells used in the 
analysis has a change in contaminant concentration specific to 
the contaminant and scenario (current concentration in micrograms 
per cubic meter as an annual average less projected concentration 
under the specified scenario in micrograms per cubic metre as an 
annual average). The change in contaminant concentration (CONC) 
for a grid cell: is multiplied by the population (PoP) ofjthevar ia 
cell. The products of these multiplications are summed over all 
grid cells. The sum’ is then multiplieavby thevcancer tna ee 
factor (URF) and divided by 70 to yield the change in annual 
Cancer mortalircy. 


Using arsenic as an example, the change in concentration 
multiplied by the population summed over all grid cells is as 
follows by scenario: 


Change in Unit Risk Change in 
. Concentration Factor Annvwal cancer 
SEAN x Population Mortality” 
Scenarios "A" & "pf" 117,335 be K LORS ey 
Scenario "B" 128,240 21/3, Ext Ox 2 79 
Scenario "c" 105,577 43. x 1072 6.5 
Scenario "E" sislice i: ak Aa Oe > PAO 


: 117 33 Seen oy TO eT. 2 


The unit risk. faéctom fom each carcinogen for which informatvor is 
available is shown im @xhibit H.1. 


H..3; NITROGEN DIOXIDE FONCTIONS 


Exposure-response functions for nitrogen dioxide have been 
developed for: | 

eye irritation; 

Visibility Tepifects 7 and 


materials damage. 


The functions Itstedt bellow are derived? rom ‘the’ Hater ial 


presented in Appendix oS 


Hoo... Bye Irmrbation 


The eye irritation damage function depends upon measurement of 
average daily high concentration of NO». We only had information 
on annual average concentration and therefore were unable to 
calculate changes eye irritation under the proposed 


implementation scenarios due to reduced emissions of NO>. 


OO  —___s 


EXHIBIT H.1: ‘. CANCER UNIT RISK FACTORS 


Contaminant Unit Risk Factor 


Acrylonitrile Sue ee LO 
Arsenic sae STOR 
Asbestos KEG rt Om 
Benzene $a ta ose 
1,3 Butadiene Pas. x 207" 
Cadmium we uo 7* 
Carbon Tetrachloride aes): xt Bone 
Chloroform 2 Bit LORS 
Chromium direc Leres 
Epichlorohydrin ae) ae oe 
Ethylene Dichloride 2.6 x 107° 
Formaldehyde phorsiner> 
Methylene Chloride AN acu 
Nickel awerse Tone 
Polychlorinated Biphenols pine savorie 
Propylene Oxide | 6.67 x 107” 
Styrene . SA7 xevos 
Toluene Diisocyanate , 6. 7S OES 
1,1,1-Trichloroethane 

(Methyl Chloroform) 1.6 x 107° 
Trichloroethylene Wien do 1b 
Vinyl Chloride . 26h Ore 


Source: Appendix C 


Mies. 2 (Vast i la ty meicecss 


We adjusted the function given in Appendix C as follows: 
. Sa, 


Total 


YY QL ie = By Xa2hSpx lini géa/tisch0n0s5 Xe NOo rer uyNOnn) ex 
Per Year LOQO/ANO5 m}s) J. eX pPOP/ HE 
Where: 
E = The mean 1986 rate of exchange (1.3894) to convert 
L986, U.oeuto. £986 Canaciian yao! lars. 
ln = The natural logarithm. 
0.035 = A mean value for the function, which varies according 
to urban: or rural: location. § (Vaking’ a mean value 


biases the result, but the difference between urban and 
rural location is not as great as the upper and lower 
BOUNGS  LOr abl hs pGUNnGEL OM t) 


NOoa = “The seurrentsNC>oeconcentration for the grid cell. 

NOop = The NO» concentration, for the same grid cell under the 
specified implementation scenario. 

POP = The population of the grid cell. 

HF = Apfactorn) (24.3)-atorconvert-population toethe inumber of 
households. | 


The function given above is the central case estimate. To get 
the lower bound estimate, the constant 205 is replaced by 100. 
In other words, the lower bound is equal to the central case 
dividedzbyo2. O05: 


As was the case for the cancer mortality equation, the 


calculation is performed for each of the over 15,000 grid cells 


Hi =*.6 


and then summed. Specifically, 


ln [1/(1+{0.035 x (NOjgq, - NOgp) x 100/NO2,})] x POP/HF 
is calculated for each grid cell. These values are then summed 
over all grid cells. And the sum is multiplied by E x 205 to get 


the final result. rast “ 


Note that the output is a dollar value; there is no physical 
value such as "number of days of poor visibility", analogous to 


the niimber of deaths, hospital days, etc? in the other functions-~ 


H.3.3 Materials Damage 


The materials damage function for nitrogen dioxide is expressed 
in similar terms. The equation describes the annual materials 
damage value per person. For the central case, the materials 


damage equation is: 


Total 


Value =F xX Esstrxrycnosn” ~ NOos) x LeO NOx) xX Por 
Per Year 


The symbols have the same meaning as in the visibility equation. 
Thus, materials damage is a function of the percentage change in 


nitrogen dioxide concentration and population. 


The function shown above is the central case. The Jowér>bound=is 
calculated as 075° times the central case: 


cS Semtaeh 


H.4 SULPHUR DIOXIDE FUNCTIONS 


H.4.1 Mortality 


An exposure-response function is available for deaths 
attributable to sulphur dioxide. The function in Appendix C is 
adjusted as follows to calculate absolute change in the number of 


annual deaths. 


Deaths = 2.672 x 107° x POP x (SOz,-SO9R) 
Where: 
SOo,n = The SO5 concentration for a grid cell under existing 
regulations. 
SOop =i—THe SOsSConcentrationt for augqrid celt undervatspecwiac 
implementation scenario. . 
POP = The population of a grid cell. 


The function presented above is the central case estimate for 
reduction in annual mortality due to reduced concentrations of 
sulphur dioxide. The low case is calculated by applying the 

lower bound estimate of the value of a statistical life to the 


calculated reduction in mortality. 


Hi4.2 sHospital Days» for Respiratory, Conditions 


The function of hospital days for respiratory conditions (HDRC) 


oS: 
HDRC = 4.0 x 107* x (SO5,.- SOQ5—5) % POP 


The variables S05,, SOop and POP have the same meaning as in the 
mortality equation. 


The lower bound estimate is obtained by replacing. 4..0 x 1622603 
the central case equation presented above with 1.5 x POFTUF In 


other words -the lower bound is: 0.375 times the "Central case. 


H.4.3 Hospital Admissions for Respiratory Disease 


The function for hospital admissions for respiratory disease 
(HARD) is: 


HARD = 4.065 x 107> x (SO, - SOop) x POP 


The variables"SO.,;f'SO5—5 and POP have the same meaning as in the 
mortality equation. 


The lower bound estimate is calculated by replacing 4.065 x 107° 
in the central case equation presented above with 0.63 x ics, 


Alternately, the lower bound is 0.155 times the central case. 
H.5 PARTICULATE MATTER FUNCTIONS 


Six exposure response functions exiac hom, OM a They relate wcas 


death; 

emergency room visits; 

restricted activity days; 

chronic obstructive pulmonary disease; 
Vile Lane 


materials damage. 


ir eo se 


The functions below are derived from those in Appendix C. In 


each case: 


PMa = Théeoriginal PM Goncentration for a griducell. 
PMp = The PM concentration for a grid cell under the 


scenario considered. 


E =» Thel méan 1986, ratevofexchange '(1.33894,) «to convert 
1936 US .tOnlOS6 .<Ganadian sdodLass, 
HE = A factor (2.8) to convert the population of a grid 


cell to the number of households. 
POP = The population of a grid cell. 


H.5.1 Death 
Deaths = 2.21 x 107© x POP x (PM,-PMg) 


This equation gives the change in annual deaths due to a change 
in particulate matter concentration. The lower bound and central 
value are calculated by applying the corresponding values of a 
statistical life to the reduction in annual mortality estimated 
using the above equation. 


H.5.2 Emergency Room Visits 
ERV = 1.3.x 1074 x POP x (PM,-PMg) 


The change in the number of emergency room visits per year, ERV, 
is related to the change in particulate matter concentration. 
The lower bound estimate of emergency room visits is calculated 
by replacing the Central case coerricient-of 1. 3> x10" by 

2 x Lorey 


H - 10 


HiS.3 Restricted Activity Davs 


RAD = 2.736 x 1072 x POP x (PMa-PMg) - Emergency Room 
Ripe ten | Visits 

The change in the number of restricted activity days, RAD, is 
reduced by the change in the number of emergency room visits 
(ERV) to avoid double counting. The lower bound estimate of 
restricted activity days is calculated by replacing the central 
case coefficient of 2.736 x 107% by 1.938 x For! +r other 
words, the lower bound is equal to 0.70833 times the central 


case. 


H.5.4 Chronic Obstructive Pulmonary Disease 


The function in Appendix C deals with changes in risk of 
pulmonary disease in adults and children. It does not predict 
actual change in, or values relating to, pulmonary disease. It 


was therefore impractical to use it in the analysis. 


H.5S35:> Visibibity Petect= 


Total 


Value i ee Oe. eu (PMa/PMp) x POP/HY 
Per Year 


The lower bound estimate is calculated as.0.5 times the central 
case. 


As with visibility effects from nitrogen dioxide, this function 


returns an annual dollar value for improved visibility. 


H.5.6 Materials Damage 


Total 


Value = Ex 7.4 x (PMa-PMp) x POP/HF 
Per Year 


jc Rear ik 


This function: also ‘returns’ an’ annual dollar value.of reduced 
materials damage. 

ca ae 
To get Pha lover bound estimate, the central case coefficient of 
7.4-is replaced’ by. 0.52. In other words, the lower bound is 


0.07027 times the central. case: 
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EXHIBIT .224.: PROFILE OF ARSENIC CONCENTRATION DUE TO REGULATION 
308 SOURCE EMISSIONS ACROSS SOUTHERN ONTARIO 
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EXHIBIT 3.16: PROFILE OF ASBESTOS CONCENTRATION DUE TO REGULATION > 
308 SOURCE EMISSIONS ACROSS SOUTHERN ONTARIO 
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EXHIBIT I.8: PROFILE OF BENZENE CONCENTRATION DUE TO REGULATION 
308 SOURCE EMISSIONS ACROSS SOUTHERN ONTARIO 
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ONTARIO | 
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PROFILE OF 2-BUTANONE CONCENTRATION DUE TO 
REGULATION 308 SOURCE EMISSIONS ACROSS SOUTHERN 
' ONTARTO 
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EXHIBIT I.14: PROFILE OF BUTYL ACRYLATE CONCENTRATION DUE TO 
REGULATION 308 SOURCE EMISSIONS ACROSS SOUTHERN 
ONTARIO 
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EXHIBIT 1.16: )\\PROFILE, OF (CADMIUM CONCENTRATION ) DUEwTOsREGULATION 
308 SOURCE EMISSIONS ACROSS SOUTHERN ONTARIO 
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EXHIBIT I.18: PROFILE OF CARBON TETRACHLORIDE CONCENTRATION, DUE 
TO REGULATION 308 SOURCE EMISSIONS ACROSS SOUTHERN 
“ONTARIO 
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EXHIBIT I.20: PROFIEE OF CDDr CONCENTRATION (DUET OIRE Gi Air Onmouo 
SOURCE EMISSIONS ACROSS SOUTHERN ONTARIO 
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EXHIBIT 1.22:  PROPTLE OF CDF CONCENTRATION. DUE TOL REGURA LON soc 
SOURCE EMISSIONS ACROSS SOUTHERN ONTARIO 
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EXHIBIT I.24: PROFILE OF CHLORINE CONCENTRATION DUE TO REGULATION 
308 SOURCE EMISSIONS ACROSS SOUTHERN ONTARIO 


Chlorine 


Concentration 
(ug/m3) 


~84 -83 -82 -81 -80 -79 -78 -77 -76 -75 -74 
Longitude 
West - East 


Health Threshold «= 30 ug/m83 


Legend 


Existing 

Scenarios A & D 

Scenario B 

Scenario C Seer ea 


Scenario E Pe atte URS 


22 


SUOISIWWS OU 


GG « 
OOF - 65 
MS fs -S 


Wy “bs “ suasuag G 


SNOISS TWA 
AOUNOS 8Oe€ NOLLVINSAY WOUT AGIXOIG ANTYOIWOD OL GHSOdXH SWVAUV 


*GC°L LIGTHXa 


LO 


8S SS 


EAU BD sec: s 


PROFILE OF CHLORINE DIOXIDE CONCENTRATION DUE TO 
REGULATION 308 SOURCE EMISSIONS ACROSS SOUTHERN 
ONTARIO 


Chlorine Dioxide 
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EXHIBIT I.28: PROFILE OF CHLOROFORM CONCENTRATION DUE TO 
REGULATION 308 SOURCE EMISSIONS ACROSS SOUTHERN 
ONTARIO 
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EXHIBIT I.30: PROFILE OF CHROMIUM CONCENTRATION DUE TO REGULATION 
308 SOURCE EMISSIONS ACROSS SOUTHERN ONTARIO 
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EXHIBIT I.32: PROFILE OF EPICHLOROHYDRIN: CONCENTRATION, DUEY TO 


REGULATION 308 SOURCE EMISSIONS ACROSS SOUTHERN 
_ ONTARIO 
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EXHIBIT 1.34: PROFILE OF ETHYL BENZENE CONCENTRATION DUE TO 
REGULATION 308 SOURCE EMISSIONS ACROSS SOUTHERN 
ONTARIO 
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EXHIBIT’I 363°“ PROFILE OF ETHYLENE’ CONCENTRATION .DUE_ TO. RECULATION 
308 SOURCE EMISSIONS ACROSS SOUTHERN ONTARIO 
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EXHIBIT I.38: PROFILE OF ETHYLENE DIBROMIDE CONCENTRATION DUE TO 


REGULATION 308 SOURCE EMISSIONS ACROSS SOUTHERN 
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EXHIBIT I.42: PROFILE OF BUTYL CELLOSOLVE CONCENTRATION DUE TO 
REGULATION 308 SOURCE EMISSIONS ACROSS SOUTHERN 
ONTARIO 
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EXHIBIT 1.44: PROFILE OF CELLOSOLVE CONCENTRATION DUE TO 
REGULATION 308 SOURCE EMISSIONS ACROSS SOUTHERN 
' ONTARIO 
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PROFILE OF CELLOSOLVE ACETATE CONCENTRATION DUE TO 
REGULATION 308 SOURCE EMISSIONS ACROSS SOUTHERN 
' ONTARIO 
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EXHIBIT, I.485 


PROFILE OF ETHYLENE OXIDE CONCENTRATION DUE TO 
REGULATION 308 SOURCE EMISSIONS ACROSS SOUTHERN 


' ONTARIO 
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EXHIBIT I.50: PROFILE OF FLUORIDES CONCENTRATION DUE TO 
REGULATION 308 SOURCE EMISSIONS ACROSS SOUTHERN 
| ONTARIO 
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EXHIBIT 1.52: \ PROFILE .OFe FORMALDEHYDE“ CONCENTRATION. DUE.ooe 
REGULATION 308 SOURCE EMISSIONS ACROSS SOUTHERN 
| ONTARIO : 
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EXHIBIT 1.54: PROFILE OF HYDROGEN CHLORIDE CONCENTRATION DUE TO 
REGULATION 308 SOURCE EMISSIONS ACROSS SOUTHERN 
ONTARIO 
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EXHIBIT I.56: PROFILE OF LEAD CONCENTRATION DUE TO REGULATION 308 
SOURCE EMISSIONS ACROSS SOUTHERN ONTARIO 
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EXHIBIT I.58: PROFILE OF MANGANESE COMPOUNDS CONCENTRATION ;DUE TO 
REGULATION 308 SOURCE EMISSIONS ACROSS SOUTHERN 
ONTARIO 
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EXHIBIT I.60: PROFILE OF MERCURY CONCENTRATION DUE TO REGULATION 
308 SOURCE EMISSIONS ACROSS SOUTHERN ONTARIO 
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EXHIBIT I.62: PROFILE OF METHYL ISOBUTYL KETONE CONCENTRATION DUE 
TO REGULATION 308 SOURCE EMISSIONS ACROSS SOUTHERN 
- ONTARIO 
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EXHIBIT 1.64: PROFILE OF METHYLENE: CHLORIDE CONCENTRAL LONGUE 1 
REGULATION 308 SOURCE EMISSIONS ACROSS SOUTHERN 
ONTARIO 
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EXHIBIT I.66: PROFILE OF NAPHTHALENE CONCENTRATION DUE TO. 
REGULATION 308 SOURCE EMISSIONS ACROSS SOUTHERN 
ONTARIO 
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EXHIBIT I.68: PROFILE OF NICKEL CONCENTRATION. DUE TO REGULATION 
308 SOURCE EMISSIONS ACROSS SOUTHERN ONTARIO 
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EXHIBIT I.70: PROFILE OF NITROGEN OXIDES CONCENTRATION DUE TO 
REGULATION 308 SOURCE EMISSIONS ACROSS SOUTHERN 
' ONTARIO 
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EXHIBIT I.72: PROFILE OF PERCHLOROETHYLENE CONCENTRATION DUE TO 
REGULATION 308 SOURCE EMISSIONS ACROSS SOUTHERN 
ONTARIO 
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EXHIBIT 1.74: PROFILE OF PHENOL CONCENTRATION DUE TO REGULATION 
308 SOURCE EMISSIONS ACROSS SOUTHERN ONTARIO 
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EXHIBIT I.76: PROFILE OF POLYCHLORINATED BIPHENYLS CONCENTRATION 

DUE TO REGULATION 308 SOURCE EMISSIONS ACROSS 
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EXHIBIT 1.78: PROFILE OF PROPTONALDERY DE CONCENTRATION SD nao 
REGULATION 308 SOURCE EMISSIONS ACROSS SOUTHERN 
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EXHIBIT I.80: PROFILE OF PROPIONIC ACID CONCENTRATION DUE TO 
REGULATION 308 SOURCE EMISSIONS ACROSS SOUTHERN 
| ONTARIO 
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EXHIBIT I.82: PROFILE OF PROPYLENE OXIDE CONCENTRATION DUE TO 
REGULATION 308 SOURCE EMISSIONS ACROSS SOUTHERN 
- ONTARIO 


weg 5 


Propylene Oxide 


Concentration 
(ug/m3) (1E-4) 


-84 -83 -82 -81 -80 -79 -78 -77 -76 -75 -74 
Longitude 
West - East 


Health Threshold = 900 ug/m3 


Legend 


Existing 

Scenarios A &D 

Scenario B LA Renee 
Scenario G@ 5. Seecee Beep 


WOOF 
— 
SUOTSIWWd OU 
QGF « 
QGF - 8S 
ii = G 
Wy “bs “¢ suasuag G 


83 


parent snes eens ae 
eteistensacsnag |) 


v 


SNOISSINA HOUNOS 8Oe€ NOTLVWIOOTA Wout WOINSTINS OL GASOdXH SVSUV +e8°L LIATHXKa 


Loses 


ES 


EXHIBIT I-84: PROFILE OF SELENIUM CONCENTRATION DUE TO REGULATION 
' 308 SOURCE EMISSIONS ACROSS SOUTHERN ONTARIO 
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EXHIBIT I.86: PROFILE OF SILICA CONCENTRATION DUE TO REGULATION 
308 SOURCE EMISSIONS ACROSS SOUTHERN ONTARIO 
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EXHIBIT 1.88: | PROFIGE. OF. STYRENE...CONCENEFRATION- DUE-TO REGULATION 
308 SOURCE EMISSIONS ACROSS SOUTHERN ONTARIO 
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EXHIBIT I.90: PROFILE OF SULPHUR DIOXIDE CONCENTRATION DUE. TO 
REGULATION 308 SOURCE EMISSIONS ACROSS SOUTHERN 
| ONTARIO 
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EXHIBIT I.92: PROFILE OF TOTAL PARTICULATE CONCENTRATION DUE TO 
REGULATION 308 SOURCE EMISSIONS ACROSS SOUTHERN 
ONTARIO 
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EXHIBIT I.94: PROFILE OF TIN CONCENTRATION DUE TO REGULATION 308 
SOURCE EMISSIONS ACROSS SOUTHERN ONTARIO 


eS ES 


oe! ” 
wr gl 


Tin 
Concentration 


(ug/m3) 
0004 = 


@) Ag) O 4 [a eseetscesneetssotmuntnceinnanesncenmaassneeincooternanntansecnt 


©) if ' ' if ‘ ” ' eat 
-84 -83 -82 -81 -80 -79 -78 -77 -76 -75 -74 
Longitude 
West - East 


Health Threshold = 2.0 ug/m3 


Legend 


Existing rn ee ees 


Scenarios A &D 


Scenario B = 
Scenario © iE hades en 
Scenario E i ae 


4 


"ROM REGULATION 308 SOURCE EMISSIONS 


5S Persons / sq. km 


50 
> 100 
no emnnisions 


pavers 


PARES ae 


ato agey 
= 


I 


« 
a 


Fe Yevee re) 
erriveg teesvaesow 
pate verteertise st xt tee 
we 
Gre 
LEE 


eee ery 
Shape spp mame bp Py 


Tecan 
Rees bee ew 
Faas yom 


EXPOSED TO TOLUENI 


Anais 
Nees cue Siena) 


AREAS 


EXHIBIT 1.95 


EXHIBIT 1.96: . PROFILE OF TOLUENE CONCENTRATION DUE TOMREGULATION 
308 SOURCE EMISSIONS ACROSS SOUTHERN ONTARIO 
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EXHIBIT I.98: PROFILE OF TOLUENE DIISOCYANATE CONCENTRATION DUE 
TO REGULATION 308 SOURCE EMISSIONS ACROSS SOUTHERN 
- ONTARIO 
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EXHIBIT 1.100: PROFILE OF ‘1,1),1-TRICHLOROETHANE (CONCENTRALOONMOUS 
TO REGULATION 308 SOURCE EMISSIONS ACROSS SOUTHERN 
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EXHIBIT I.102: PROFILE OF TRICHLOROBENZENES CONCENTRATION DUE TO 
REGULATION 308 SOURCE EMISSIONS ACROSS SOUTHERN 
| ONTARIO 
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EXHIBIT I.104: PROFILE OF TRICHLOROETHYLENE CONCENTRATION DEETIO 
REGULATION 308 SOURCE EMISSIONS ACROSS SOUTHERN 
ONTARIO 
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EXHIBIT I.106: PROFILE OF VINYL CHLORIDE CONCENTRATION DUE TO 
REGULATION 308 SOURCE EMISSIONS ACROSS SOUTHERN 
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EXHIBIT I.108: PROFILE OF XYLENES CONCENTRATION DUE TO REGULATION 
308 SOURCE EMISSIONS ACROSS SOUTHERN ONTARIO 
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